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official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


1  Introduction 


Estimating  nearshore  wind-wave  growth  and  transformation  is  a  critical 
component  of  most  coastal  engineering  projects,  e.g.,  predicting  bathymetric  and 
shoreline  change,  estimating  navigation  channel  shoaling  and  migration, 
designing  or  repairing  coastal  structures,  assessing  navigation  conditions,  and 
evaluating  natural  evolution  of  coastal  inlets  or  beaches  versus  consequences  of 
engineering  actions.  Nearshore  wave  propagation  is  influenced  by  complex 
bathymetry  (including  shoals  and  navigation  channels);  tide-,  wind-,  and  wave¬ 
generated  currents;  tide-  and  surge-induced  water  level  variation;  and  coastal 
structures.  Use  of  numerical  wave  models  has  become  widespread  to  represent 
wave  transformation  primarily  because  of  their  increasing  sophistication  and 
economy  of  application  relative  to  the  large  expense  of  field  measurements  or 
physical  model  studies. 

This  report  describes  version  3.0  of  the  steady-state  spectral  wave  model, 
STWAVE.  The  purpose  of  STWAVE  is  to  provide  an  easy-to-apply,  flexible, 
and  robust  model  for  nearshore  wind-wave  growth  and  propagation.  Upgrades  to 
the  model  since  the  release  of  version  2.0  (Smith,  Resio,  and  Zundel  1999) 
include: 

a.  Calculation  of  radiation  stress  gradients.  These  stresses  can  be  used  to 
calculate  wave-driven  currents  and  water  level  changes  in  a  circulation 
model.  This  feature  has  been  used  to  link  with  the  Advanced  3-D 
Circulation  model  (ADCIRC)  (e.g.,  Militello  et  al.  2000;  Smith  and 
Ebersole  2000) 

b.  Increased  model  efficiency  of  25  to  40  percent  compared  to  STWAVE 
version  2.0.  The  greatest  speedup  is  for  cases  that  include  wave-current 
interaction. 

c.  Addition  of  an  input  current  option  ICUR  =  2  to  read  a  single  current  field 
for  multiple  input  spectra. 

d.  Addition  of  an  output  file  designating  grid  cells  where  wave  breaking 
occurs. 

e.  Revised  format  of  the  spectral  output  file.  This  modification  eliminates 
an  error  in  free  format  reads  for  spectra  with  high  energy  levels. 

f.  Revised  calculation  of  wave  parameters  on  the  offshore  boundary  that 
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occurred  in  version  2.0  when  the  offshore  boundary  included  both  land 
and  water  cells. 

g.  Improved  Surface-Water  Modeling  System  (SMS)  user  interface  and 
supporting  software  for  grid  generation,  interpolation  of  current  fields, 
generation  of  input  spectra,  and  visualization  of  model  output. 

Version  7.1  of  SMS  provides  significant  improvements  for  conversion  of 
coordinate  systems,  interpolating  current  fields,  and  generating  and 
visualizing  input/output  wave  spectra.  This  version  also  includes  the 
additional  input  parameters  required  for  STWAVE  version  3.0  and 
visualization  of  breaking  wave  fields. 

This  report  describes  procedures  for  using  STWAVE.  In  Chapter  2,  an 
overview  of  the  model  governing  equations  and  the  numerical  discretization  is 
presented.  In  Chapter  3,  the  necessary  input  files  are  described.  Guidelines  for 
selecting  parameter  values  are  also  given.  In  Chapter  4,  the  content  and  format 
of  the  STWAVE  output  files  are  described.  Chapter  5  gives  a  tutorial  on  SMS 
application  (version  7.1)  using  data  from  the  Yaquina  Entrance  in  Oregon. 
Chapter  6  provides  three  example  applications  of  the  model  for  wave  growth  and 
transformation  over  simple  bathymetry  and  wave  transformation  over  complex 
bathymetry  with  wave-current  interaction.  Sample  input  and  output  files  are 
given  in  Appendices  A-G. 
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2  Governing  Equations  and 
Numerical  Discretization 


This  chapter  gives  an  overview  of  the  phase-averaged  spectral  wave  model 
STWAVE  (STeady-state  spectral  WAVE  model)  (Resio  1987, 1988a,  1988b; 
Davis  1992;  Smith,  Resio,  and  Zundel  1999).  STWAVE  is  a  steady-state  finite- 
difference  model  based  on  the  wave  action  balance  equation.  This  report 
describes  STWAVE  version  3.0,  which  includes  calculation  of  radiation  stresses 
and  identification  of  regions  of  active  wave  breaking. 


Model  Capabilities 

The  purpose  of  applying  nearshore  wave  transformation  models  is  to 
describe  quantitatively  the  change  in  wave  parameters  (wave  height,  period, 
direction,  and  spectral  shape)  between  the  offshore  and  the  nearshore  (typically 
depths  of  40  m  or  less).  In  relatively  deep  water,  the  wave  field  is  fairly 
homogeneous  on  the  scale  of  kilometers;  but  in  the  nearshore,  where  waves  are 
strongly  influenced  by  variations  in  bathymetiy,  water  level,  and  current,  wave 
parameters  may  vary  significantly  on  the  scale  of  tens  of  meters.  Offshore  wave 
information  is  typically  available  from  a  wave  gauge  or  a  global-  or  regional- 
scale  wave  hindcast  or  forecast.  Nearshore  wave  information  is  required  for  the 
design  of  almost  all  coastal  engineering  projects.  Waves  drive  sediment  transport 
and  nearshore  currents,  induce  wave  setup  and  runup,  excite  harbor  oscillations, 
or  impact  coastal  structures.  The  longshore  and  cross-shore  gradients  in  wave 
height  and  direction  can  be  as  important  as  the  magnitude  of  these  parameters  for 
some  coastal  design  problems. 

STWAVE  simulates  depth-induced  wave  refraction  and  shoaling,  current- 
induced  refraction  and  shoaling,  depth-  and  steepness-induced  wave  breaking, 
diffraction,  wind-wave  growth,  and  wave-wave  interaction  and  whitecapping  that 
redistribute  and  dissipate  energy  in  a  growing  wave  field. 

A  wave  spectrum  is  a  statistical  representation  of  a  wave  field. 

Conceptually,  a  spectrum  is  a  linear  superposition  of  monochromatic  waves.  A 
spectrum  describes  the  distribution  of  wave  energy  as  a  function  of  frequency 
(one-dimensional  spectrum)  or  frequency  and  direction  (two-dimensional 
spectrum).  An  example  of  a  one-dimensional  wave  spectrum  is  given  in 
Figure  1 .  The  peak  period  of  the  spectrum  is  the  reciprocal  of  the  frequency  of 
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Figure  1 .  Sample  one-dimensional  wave  spectrum 

the  peak  of  the  spectrum.  The  wave  height  (significant  or  zero-moment  wave 
height)  is  equal  to  four  times  the  square  root  of  the  area  under  the  spectrum.  For 
the  example  spectrum  given  in  Figure  1,  the  peak  frequency  is  0.105  Hz,  the  peak 
period  is  9.5  sec,  and  the  wave  height  is  2.8  m.  STWAVE  is  based  on  the 
assumption  that  the  relative  phases  of  the  spectral  components  are  random,  and 
thus  phase  information  is  not  tracked  (i.e.,  it  is  a  phase-averaged  model).  In 
practical  applications,  wave  phase  information  throughout  a  model  domain  is 
rarely  known  accurately  enough  to  initiate  a  phase-resolving  model.  Typically, 
wave  phase  information  is  only  required  to  resolve  wave-height  variations  near 
coastal  structures  for  detailed,  near-field  reflection  and  diffraction  patterns. 

Thus,  for  these  situations,  a  phase-resolving  model  should  be  applied. 


Model  Assumptions 

The  assumptions  made  in  STWAVE  version  3.0  are: 

a.  Mild  bottom  slope  and  negligible  wave  reflection.  STWAVE  is  a  half¬ 
plane  model,  meaning  that  wave  energy  can  propagate  only  from  the 
offshore  toward  the  nearshore  (+87.5  deg  from  the  x-axis  of  the  grid, 
which  is  typically  the  approximate  shore-normal  direction).  Waves 
reflected  from  the  shoreline  or  from  steep  bottom  features  travel  in 
directions  outside  this  half  plane  and  thus  are  neglected. 
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Forward-scattered  waves,  e.g.,  waves  reflected  off  a  structure  but 
traveling  in  the  +x-direction,  are  also  neglected. 

b.  Spatially  homogeneous  offshore  wave  conditions.  The  variation  in  the 
wave  spectrum  along  the  offshore  boundary  of  a  modeling  domain  is 
rarely  known,  and  for  domains  on  the  order  of  tens  of  kilometers,  is 
expected  to  be  small.  Thus,  the  input  spectrum  in  STWAVE  is  constant 
along  the  offshore  boundary.  Future  versions  of  the  model  will  allow 
variable  input. 

c.  Steady-state  waves,  currents,  and  winds.  STWAVE  is  formulated  as  a 
steady-state  model.  A  steady-state  formulation  reduces  computation  time 
and  is  appropriate  for  wave  conditions  that  vary  more  slowly  than  the 
time  it  takes  for  waves  to  transit  the  computational  grid.  For  wave 
generation,  the  steady-state  assumption  means  that  the  winds  have 
remained  steady  sufficiently  long  for  the  waves  to  attain  fetch-limited  or 
fully  developed  conditions  (waves  are  not  limited  by  the  duration  of  the 
winds). 

d.  Linear  refraction  and  shoaling.  STWAVE  incorporates  only  linear  wave 
refraction  and  shoaling,  thus  does  not  represent  wave  asymmetry.  Model 
accuracy  is  therefore  reduced  (wave  heights  are  underestimated)  at  large 
Ursell  numbers. 

e.  Depth-uniform  current.  The  wave-current  interaction  in  the  model  is 
based  on  a  current  that  is  constant  through  the  water  column.  If  strong 
vertical  gradients  in  current  occur,  their  modification  of  refraction  and 
shoaling  is  not  represented  in  the  model.  For  most  applications,  three- 
dimensional  current  fields  are  not  available. 

f  Bottom  friction  is  neglected.  The  significance  of  bottom  friction  on  wave 
dissipation  has  been  a  topic  of  debate  in  wave  modeling  literature. 

Bottom  friction  has  often  been  applied  as  a  tuning  coefficient  to  bring 
model  results  into  alignment  with  measurements.  Although  bottom 
friction  is  easy  to  apply  in  a  wave  model,  determining  the  proper  friction 
coefficients  is  difficult.  Also,  propagation  distances  in  a  nearshore  model 
are  relatively  short  (tens  of  kilometers),  so  that  the  cumulative  bottom 
friction  dissipation  is  small.  For  these  reasons,  bottom  friction  is 
neglected  in  STWAVE. 

g.  Linear  radiation  stress.  Radiation  stress  is  calculated  based  on  linear 
wave  theory. 

Ongoing  research  will  enhance  present  model  capabilities  and  eliminate  some 
model  assumptions.  The  following  sections  describe  wave  propagation  and 
source/sink  terms  in  STWAVE  version  3.0. 
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Governing  Equations 

Interaction  of  waves  with  currents  is  considered  in  a  reference  frame  moving 
with  the  current.  Wave  parameters  in  this  frame  are  denoted  with  the  subscript  r, 
for  being  “relative”  to  the  current,  and  parameters  in  the  nonmoving  reference 
frame  are  subscripted  a,  for  “absolute.”  The  wave  dispersion  relationship  is 
given  in  the  moving  reference  frame  as  (Jonsson  1990  and  others): 

col  =  gk  tanh  kd  (1) 

where 

co  =  angular  frequency 
g  =  gravitational  acceleration 
k  =  wave  number 
d  =  water  depth 

In  the  absolute  frame  of  reference,  the  dispersion  equation  is: 

coa=cor  +  kU cos (5  -a)  (2) 

where 

U  =  current  speed 

S  =  direction  of  the  current  relative  to  a  reference  frame  (the  x-axis,  here) 
a  =  wave  orthogonal  direction  (normal  to  the  wave  crest)  (Figure  2) 


The  wave  number  is  solved  by  substituting  Equation  1  into  Equation  2  and 
iteratively  solving  for  k.  The  wave  number  and  wavelength  (L=(2n)lk)  are  the 
same  in  both  reference  frames. 

Solutions  for  refraction  and  shoaling  also  require  wave  celerities,  C,  and 
group  celerities,  Cg,  in  both  reference  frames.  In  the  reference  frame  relative  to 
the  current, 


c. 


k 


C. 


Sr 


0.5  C. 


1  + 


V 


2  kd  ^ 
sinh2M  j 


(3) 

(4) 
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Figure  2.  Definition  sketch  of  wave  and  current  vectors 


The  direction  of  both  the  relative  celerity  and  group  celerity  is  a,  the  wave 
orthogonal  direction.  In  the  absolute  reference  frame, 


ca  =  Cr+Ucos(S  -a) 

(5) 

(6) 

where  subscript  i  is  tensor  notation  for  the  x-  and  y-components.  The  direction  of 
the  absolute  celerity  is  also  in  the  wave  orthogonal  direction.  The  absolute  group 
celerity  defines  the  direction  of  the  wave  ray,  so  the  wave  ray  direction  (Figure  2) 
is  defined  as: 


=  tan 


f  Cgr  sina  +  U sin<5 
Cgrcosa  +  U  cosS 


\ 

J 


(7) 


The  distinction  between  the  wave  orthogonal  (direction  perpendicular  to  the 
wave  crest)  and  the  wave  ray  (direction  of  energy  propagation)  is  important  in 
describing  wave-current  interaction.  Without  currents,  the  wave  rays  and 
orthogonals  are  the  same,  but  with  currents,  the  wave  energy  moves  along  the 
rays  whereas  the  wave  direction  is  defined  by  the  orthogonals. 

The  wave  orthogonal  direction  for  steady-state  conditions  is  given  by  (Mei 
1989;  Jonsson  1990): 


,  Da  Crk  Dd_  k,  DU, 
80  DR  sinh  2kd  Dn  k  Dn 
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where 


D  =  derivative 

R  =  coordinate  in  the  direction  of  the  wave  ray 
n  =  coordinate  normal  to  the  wave  orthogonal 

The  governing  equation  for  steady-state  conservation  of  spectral  wave  action 
along  a  wave  ray  is  given  by  (Jonsson  1990): 

N  8  CaCgacos(/*-a)£(roa,a)_v  S  ,a> 

V^ga/i  ^  ^  ' 

OXl  (Or  (Or 

where 

E  =  wave  energy  density  divided  by  (pw  g),  where  pw  is  density  of 
water 

S  =  energy  source  and  sink  terms 


Refraction  and  shoaling 

Refraction  and  shoaling  are  implemented  in  STWAVE  by  applying  the 
conservation  of  wave  action  along  backward  traced  wave  rays.  Rays  are  traced 
in  a  piecewise  manner,  from  one  grid  column  to  the  next.  The  two-dimensional 
wave  spectra  are  set  as  input  along  the  first  grid  column  (the  offshore  boundary). 
For  a  point  on  the  second  grid  column,  the  spectrum  is  calculated  by  back  tracing 
a  ray  for  each  frequency  and  direction  component  of  the  spectrum.  The  ray 
direction,  //,  is  determined  by  Equation  7.  Only  ray  directions  propagating 
toward  the  shore  (-87.5  to  +87.5  deg)  are  included.  Energy  propagating  toward 
the  offshore  is  neglected. 

The  wave  ray  is  traced  back  to  the  previous  grid  column,  and  the  length  of 
the  ray  segment  DR  is  calculated.  Derivatives  of  depth  and  current  components 
normal  to  the  wave  orthogonal  are  estimated  (based  on  the  orthogonal  direction 
at  column  2)  and  substituted  into  Equation  8  to  calculate  the  wave  orthogonal 
direction  at  column  1 .  Then,  the  wave  number,  wave  and  group  celerities,  and 
ray  angle  in  the  previous  column  are  calculated.  The  energy  is  calculated  as  a 
weighted  average  of  energy  between  the  two  adjacent  grid  points  in  the  column 
and  the  direction  bins.  The  energy  density  is  corrected  by  a  factor  that  is  the  ratio 
of  the  5-deg  standard  angle  band  width  to  the  width  of  the  back-traced  band  to 
account  for  the  different  angle  increment  in  the  back-traced  ray.  The  shoaled  and 
refracted  wave  energy  in  column  2  is  then  calculated  from  the  conservation  of 
wave  action  along  a  ray  (Equation  9). 

In  a  strong  opposing  current  (e.g.,  ebb  currents  at  an  entrance),  waves  may  be 
blocked  by  the  current.  Blocking  occurs  if  there  is  no  solution  to  the  dispersion 
equation  (Equation  2).  Or,  to  state  it  another  way,  blocking  occurs  if  the  relative 
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wave  group  celerity  is  smaller  than  the  magnitude  of  the  opposing  current,  so 
wave  energy  cannot  propagate  against  the  current.  In  deep  water,  blocking 
occurs  for  an  opposing  current  with  magnitude  greater  than  one-fourth  the 
deepwater  wave  celerity  without  current  (0.25  g  TJQ.it),  where  Ta  is  the  absolute 
wave  period).  If  blocking  occurs,  the  wave  energy  is  dissipated  through  breaking. 
Lai,  Long,  and  Huang  (1989)  performed  laboratory  experiments  that  showed  that 
wave  energy  can  pass  through  the  linear  blocking  point  through  nonlinear  energy 
transfers  to  lower  frequencies  (which  are  not  blocked).  These  nonlinear  energy 
transfers  are  not  included  in  STWAVE. 


Diffraction 

Diffraction  is  included  in  STWAVE  in  a  simple  manner  through  smoothing 
of  wave  energy.  The  model  smoothes  energy  in  a  given  frequency  and  direction 
band  using  the  following  form: 

Ej  {coa  ,a)  =  0.55  Ej  ( coa,a)  +  0.225[EJ+l  (coa  ,a )  +  Ej_ ,  (©„,«)]  (10) 

where  E  is  the  energy  density  in  a  given  frequency  and  direction  band,  and  the 
subscript  j  indicates  the  grid  row  index  (alongshore  position).  Equation  10 
provides  smoothing  of  strong  gradients  in  wave  height  that  occur  in  sheltered 
regions,  but  provides  no  turning  of  the  waves.  This  formulation  is  grid-spacing 
dependent,  which  is  a  serious  weakness.  Efforts  are  ongoing  to  implement  a 
more  rigorous  diffraction  representation. 


Source/sink  terms 

Surf-zone  wave  breaking.  The  wave-breaking  criterion  applied  in  the  first 
version  of  STWAVE  was  a  function  of  the  ratio  of  wave  height  to  water  depth: 

^^=0.64  (11) 

d 

where  Hmo  is  the  energy-based  zero-moment  wave  height.  At  a  coastal  entrance, 
where  waves  steepen  because  of  the  wave-current  interaction,  wave  breaking  is 
enhanced  because  of  the  increased  steepening.  Smith,  Resio,  and  Vincent 
(1997)  performed  laboratory  measurements  of  irregular  wave  breaking  on  ebb 
currents  and  found  that  a  breaking  relationship  in  the  form  of  the  Miche  criterion 
(1951)  was  simple,  robust,  and  accurate: 

Hmn  =  O.lLtanhkd  (12) 

(see  also  Battjes  1982  and  Battjes  and  Janssen  1978).  Equation  12  is  applied  in 
version  3.0  of  STWAVE  as  a  maximum  limit  on  the  zero-moment  wave  height. 
The  energy  in  the  spectrum  is  reduced  at  each  frequency  and  direction  in 
proportion  to  the  amount  of  pre-breaking  energy  in  each  frequency  and  direction 
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band.  Nonlinear  transfers  of  energy  to  high  frequencies  that  occur  during 
breaking  are  not  represented  in  the  model.  Model  grid  cells  where  the  wave 
height  is  limited  by  Equation  12  are  flagged  as  actively  breaking  cells.  These 
breaking  regions  can  be  visualized  in  SMS. 

Wind  input.  Waves  grow  through  the  transfer  of  momentum  from  the  wind 
field  to  the  wave  field.  The  flux  of  energy,  Fim  into  the  wave  field  in  STWAVE 
is  given  by  (Resio  1988a): 

Fin=x£°-0.S5Cm—  (13) 

Pw  S 

where 

X  -  partitioning  coefficient  that  represents  the  percentage  of  total 
atmosphere  to  water  momentum  transfer  that  goes  directly 
into  the  wave  field  (0.75) 

pa  -  density  of  air 

Cm  =  mean  wave  celerity 

u*  =  friction  velocity  (equal  to  the  product  of  the  wind  speed,  U,  and 
the  square  root  of  the  drag  coefficient,  CD  =  .0012+.000025LO 

In  deep  water,  STWAVE  provides  a  total  energy  growth  rate  that  is  consistent 
with  Hasselmann  et  al.  (1973). 

The  energy  gain  to  the  spectrum  is  calculated  by  multiplying  the  energy  flux 
by  the  equivalent  time  for  the  wave  to  travel  across  a  grid  cell: 


A  t  = 


Ax 

pCg  cos  am 


(14) 


where 


At  =  equivalent  travel  time 
Ax  =  grid  spacing 

(3  =  factor  equal  to  0.9  for  wind  seas 
Cg  =  average  group  celerity  of  the  spectrum 
am  =  mean  wave  direction,  relative  to  the  grid 
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Because  STWAVE  is  a  half-plane  model,  only  winds  blowing  toward  the  shore 
(+*  direction)  are  included.  Wave  damping  by  offshore  winds  and  growth  of 
offshore-traveling  waves  are  neglected. 

Wave-wave  interaction  and  whitecapping.  As  energy  is  fed  into  the  waves 
from  the  wind,  it  is  redistributed  through  nonlinear  wave-wave  interaction. 
Energy  is  transferred  from  the  peak  of  the  spectrum  to  lower  frequencies 
(decreasing  the  peak  frequency  or  increasing  the  peak  period)  and  to  high 
frequencies  (where  it  is  dissipated). 

In  STWAVE,  the  frequency  of  the  spectral  peak  is  allowed  to  increase  with 
fetch  (or  equivalently  propagation  time  across  a  fetch).  The  equation  for  this  rate 
of  change  of  fp  is  given  by: 


(fP)M 


9 

Jg 


(  V/3 

'  u.  ' 


(15) 


where  the  i  and  i+1  subscripts  refer  to  the  grid  column  indices  within  STWAVE 
and  £  is  a  dimensionless  constant  (Resio  and  Perrie  1989).  The  energy  gained  by 
the  spectrum  is  distributed  within  frequencies  on  the  forward  face  of  the 
spectrum  (frequencies  lower  than  the  peak  frequency)  in  a  manner  that  retains  the 
self-similar  shape  of  the  spectrum. 

Wave  energy  is  dissipated  (most  notably  in  an  actively  growing  wave  field) 
through  energy  transferred  to  high  frequencies  and  dissipated  through  wave 
breaking  (whitecapping)  and  turbulent/viscous  effects.  There  is  a  dynamic 
balance  between  energy  entering  the  wave  field  because  of  wind  input  and 
energy  leaving  the  wave  field  because  of  nonlinear  fluxes  to  higher  frequencies 
(Resio  1987,  1988a).  The  energy  flux  to  high  frequencies  is  represented  in 
STWAVE  as: 


rv  = 


1/2  r>3  .9/2 

gg  E,o,kp 

tanh3/4(£p<7) 


(16) 


(Resio  1987),  where 
rE  =  energy  flux 
s  =  coefficient  equal  to  30 

E,„,  total  energy  in  the  spectrum  divided  by  (pw  g) 

kp  =  wave  number  associated  with  the  peak  of  the  spectrum 

The  energy  loss  from  the  spectrum  is  calculated  by  multiplying  the  energy 
flux  by  the  equivalent  time  for  the  wave  to  travel  across  a  grid  cell  (At) 
(Equation  14)  with  /?  equal  to  1 .0  for  the  swell  portion  of  the  spectrum  and 
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0.9  for  the  sea  portion  of  the  spectrum.  This  dissipation  is  only  applied  in  the 
model  if  wind  input  is  included. 


Radiation  stress  gradients 

Gradients  in  radiation  stress  are  calculated  in  STWAVE  to  provide  wave 
forcing  to  external  circulation  models  to  drive  nearshore  currents  and  water  level 
changes  (i.e..  wave  setup  and  setdown).  Wave-driven  currents  are  generally  the 
dominant  forcing  for  sediment  transport  in  the  surf  zone.  Radiation  stress  tensors 
are  calculated  based  on  linear  wave  theory: 


s-  -  ^  if£(/'  “M1 + S(cos2a + 1)-0-5 


dfda 


(17) 


$xy  PwS  jj" 


0.5 


1  +  - 


2  kd 


sinh2fo/ 


sin  2a 


dfda 


(18) 


'yy 


P*g 


0.5 


1+- 


2  kd 


sinh  2kd 


(sin2a  +  l)-0.5 


dfda 


The  gradients  in  radiation  stress  are  calculated  as: 


(19) 


as „  a?, 

dx  dy 


(20) 


T 


y 


dx  dy 


(21) 


Values  of  zx/pw  and  zfpw  are  output  from  STWAVE  for  use  in  circulation 
modeling. 


Numerical  Discretization 

STWAVE  is  a  finite-difference  numerical  model,  formulated  on  a  Cartesian 
grid.  Grid  cells  are  square  (Ac  =  Ay).  Variable  grid  resolution  can  be  obtained 
by  nesting  model  runs.  This  is  accomplished  by  running  the  model  at  a  coarse 
resolution  and  saving  a  spectrum  at  a  nearshore  point.  This  nearshore  spectrum 
can  then  be  used  as  a  boundary  condition  for  another  .grid  of  finer  resolution.  A 
schematic  of  a  grid  is  shown  in  Figure  3.  STWAVE  operates  in  a  local 
coordinate  system,  with  the  x-axis  oriented  in  the  cross-shore  direction  (origin 
offshore)  and  the  y-axis  oriented  alongshore,  forming  a  right-handed  coordinate 
system.  The  orientation  of  the  x-axis  (±87.5  deg)  defines  the  half  plane  that  is 
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represented  in  the  model.  The  y-axis  is  typically  aligned  with  the  bottom 
contours.  Wave  angles  are  defined  in  a  mathematical  sense,  measured 
counterclockwise  from  the  x-axis. 

Lateral  boundaries  in  the  model  can  be  specified  as  land  or  water  by 
specifying  the  cell  depths  as  positive  (water)  or  negative  (land).  Note  that 
specifying  land  around  the  entire  grid  will  give  different  results  than  if  the  lateral 
boundaries  are  water.  Land  boundaries  reduce  wave  growth  near  the  boundary 
because  they  “block”  propagation  from  landward  directions.  If  the  boundaries 
are  specified  as  water,  a  zero-gradient  type  of  boundary  is  applied  that  allows 
energy,  consistent  with  neighboring  cells,  to  propagate  into  or  out  of  the  domain 
along  the  lateral  boundary. 
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3  Input  File  Description 


STWAVE  has  four  input  files.  These  files  specify  model  parameters, 
bathymetry,  incident  wave  spectra,  and  current  fields.  The  current  field  file  is 
optional  and  required  only  if  the  wave-current  interaction  is  specified.  The  other 
three  files  are  required.  The  input  files  are  most  easily  generated  using  the  SMS 
interface  discussed  in  Chapter  5,  but  may  also  be  generated  using  an  editor  or 
other  programs.  File  formats  for  the  input  files  must  remain  the  same,  regardless 
of  how  the  files  are  generated.  Parameter  names  starting  with  “N,”  “I,”  or  “J” 
indicate  integer  parameter  values  are  required.  All  other  names  indicate  real 
parameter  values.  A  schematic  of  the  STWAVE  input  and  output  files  is  given  in 
Figure  4. 


Model  Input 


Model  Parameters 


Bathymetry 


> 

> 


Incident  Wave  Spectra 
Wind,  and  Water  Levels 


Current  Fields 
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W 

A 
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Model  Output 


Selected  Wave 


U  Selected  Wave 


Fields  of  Wave  Height, 
Period,  and  Direction 


Fields  of  Radiation 
Stress  Gradients 


Fields  of  Breaker  Indices 


Figure  4.  STWAVE  input  and  output  files 

Model  Parameters 

The  model -parameter  file  specifies  options  for  running  STWAVE  and  special 
output  points.  The  first  line  of  the  model  parameter  file,  read  using  free  format 
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(parameters  must  be  delimited  by  spaces,  commas,  or  new  lines),  includes  (in  the 
following  order): 

IPRP  =  Switch  for  including  propagation  only  (IPRP  =  1)  or  both 
propagation  and  source  terms  (IPRP  =  0).  Surf-zone  wave  breaking  is  included 
for  either  option,  but  wind-wave  generation,  wave-wave  interactions,  and 
whitecapping  are  included  only  for  IPRP  =  0.  In  applications  where  the  wave 
propagation  distances  are  short  (on  the  order  of  kilometers),  propagation  only  is 
usually  sufficient.  For  longer  wave  propagation  distances  or  cases  where  only 
locally  generated  waves  exist,  the  source  terms  should  be  included. 

ICUR  =  Switch  for  including  (ICUR  =  1)  or  excluding  (ICUR  =  0)  wave- 
current  interaction.  A  current  field  must  be  specified  for  each  incident  wave 
condition  (see  “Current  Fields”  in  a  following  section)  if  ICUR  =  1 .  A  new 
option  in  STWAVE  version  3  is  ICUR  =  2  to  specify  a  single  current  field  for  all 
incident  wave  conditions.  This  option  is  useful  for  climatic  studies  where  peak 
ebb  or  peak  flood  current  fields  may  be  used  with  a  number  of  typical  or  extreme 
incident  wave  conditions.  Currents  can  increase  (opposing  current)  or  decrease 
(following  current)  wave  height  as  well  as  change  wave  direction  because  of 
wave-current  interaction.  These  interactions  are  strongest  for  short  wave  periods 
and  large  current  magnitudes.  Smith,  Militello,  and  Smith  (1998)  showed  that 
current  magnitudes  less  than  about  1  m/sec  had  minor  impacts  on  wave 
transformation  at  an  East  Coast  entrance  in  Florida  (peak  wave  periods  in  the 
range  of  5-15  sec).  Currents  are  read  from  the  file  current .  in  (or 
pro  j  ect .  cur  in  SMS). 

IBREAK  =  Switch  for  printing  (IBREAK  =  1)  or  not  printing  (IBREAK  =  0) 
binary  indices  to  indicate  cells  in  which  the  waves  are  actively  dissipating  wave 
energy.  These  breaker  indices  are  written  to  a  file  called  break 
(project .  brk  in  SMS).  Wave  breaking  is  defined  by  Equation  12.  The  value 
of  IBREAK  has  no  impact  on  the  output  wave  fields,  but  determines  if  the 
break  output  file  is  written. 

IRS  =  Switch  for  calculating  (IRS  =  1)  or  not  calculating  (IRS  =  0)  radiation 
stress  gradients.  The  radiation  stress  gradients  are  written  to  the  file 
radstress.  The  radiation  stress  gradients  are  defined  by  Equations  17-21. 

The  radiation  stress  output  file  is  not  used  by  STWAVE,  but  it  can  be  used  as 
input  to  a  circulation  model  (e.g.,  ADCIRC)  to  force  wave-driven  currents  and 
water  level  changes.  The  output  file  can  be  relatively  large  (tv  and  rv  are  written 
for  each  grid  cell  for  each  input  wave  condition). 

NSELCT  =  Number  of  special  output  points.  The  two-dimensional  wave 
spectra  will  be  saved  only  for  these  output  points  (spec .  out).  Also,  a  file  that 
summarizes  wave  height,  period,  and  direction  for  these  points  will  be  written 
(“selected  heights”  file  selhts).  Wave  height,  period,  and  direction  are  saved 
by  the  model  at  all  grid  points  (see  Chapter  4),  but  this  summary  file  makes  it 
easy  to  spot-check  results  at  critical  locations  or  to  isolate  only  the  locations  of 
interest.  The  spectral  output  file  can  become  large  if  NSELCT  is  large  and  many 
input  wave  conditions  are  run. 
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The  final  information  in  the  model-parameter  file  is  specification  of  the 
special  output  points.  These  points  are  I,  J  (x  grid  cell  index,  y  grid  cell  index) 
pairings  of  these  special  output  points.  There  must  be  NSELCT  pairs  of 
numbers.  The  values  are  again  read  in  free  format  and  may  occupy  as  many  lines 
as  needed  (one  pair  per  line). 

The  default  name  of  the  model  parameter  file  is  options .  std.  In  SMS, 
the  file  is  specified  with  a  project  name  and  the  “std”  extension,  project .  std.  A 
sample  model  parameter  file  is  given  in  Appendix  A.  The  sample  file  specifies 
that  the  model  run  will  consider  propagation  and  source  terms  (IPRP  =  0),  wave- 
current  interaction  is  included  (ICUR  =  1),  breaker  indices  are  written  to  a  file 
(IBREAK=1),  and  radiation  stress  gradients  are  written  to  a  file  (IRS=1).  Three 
output  points  have  been  selected  (at  I,  J  coordinates  (110,157),  (131,117),  and 
(146,86)). 


Bathymetry 

The  bathymetry  file  gives  the  STWAVE  grid  size  and  grid  spacing  and 
defines  the  grid  bathymetry.  The  first  line  of  the  file,  read  using  free  format, 
includes  (in  the  following  order): 

NI  =  Number  of  cross-shore  (x)  grid  cells  or  columns.  The  value  of  NI, 
together  with  the  grid  spacing  (DXINC),  determines  the  cross-shore  extent  of  the 
modeling  domain  and  the  location  of  the  offshore  grid  boundary.  The  offshore 
boundary  is  typically  placed  in  a  region  where  the  bathymetry  contours  are 
relatively  straight  and  parallel  and  where  incident  wave  information  is  available 
from  a  gauge  or  hindcast. 

NJ  =  Number  of  alongshore  (y)  grid  cells  or  rows.  The  value  of  NJ,  together 
with  the  grid  spacing,  determines  the  alongshore  extent  of  the  modeling  domain 
and  the  location  of  the  lateral  grid  boundaries.  The  lateral  extent  of  the  grid  is  a 
balance  between  keeping  the  lateral  boundaries  as  far  away  as  possible  from  the 
region  of  engineering  interest  and  keeping  the  grid  domain  small  for 
computational  efficiency.  The  grid  should  extend  sufficiently  far  alongshore  to 
include  the  bathymetric  features  that  influence  the  engineering  project  (e.g.,  if  a 
project  is  centered  on  a  coastal  entrance,  the  entire  ebb  shoal  and  any  offshore 
features  that  influence  refraction  should  be  included  in  the  model  domain). 

DXINC  =  Grid  spacing  in  meters,  which  is  the  same  in  the  x  and  y  directions. 
The  grid  spacing  determines  how  finely  bathymetric  features  and  current  fields 
are  resolved  in  the  model  and  the  spatial  resolution  of  the  output  fields  of  wave 
height,  period,  and  direction.  Grid  spacing  should  be  sufficient  to  define 
bathymetric  features  and  gradients  (shoals,  canyons,  and  channels)  that  are 
important  to  the  engineering  project.  Typical  resolutions  are  tens  to  hundreds  of 
meters.  For  a  fixed  grid  domain,  increasing  the  grid  resolution  increases 
computational  time. 

The  remainder  of  the  file  gives  the  water  depth  (relative  to  some  datum,  such 
as  mean  sea  level  or  mean  lower  low  water)  for  each  grid  cell  in  an  STWAVE 
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application,  in  meters.  Water  depths  are  defined  as  positive  numbers  and  land 
elevations  are  negative  numbers.  The  depths  are  read  with  a  free  format.  The 
file  begins  with  the  depth  at  cell  (1 ,  NJ)  and  reads  in  the  cross-shore  direction  (I 
=  1  to  NI).  The  read  is  repeated  for  J  =  NJ-1  (I  =  1  to  NI),  and  progresses  to  J  = 
1.  The  default  name  for  the  bathymetry  file  is  depth,  in.  In  SMS,  the  file  is 
specified  as  pro  j  ect .  dep.  A  sample  bathymetry  file  is  given  in  Appendix  B 
The  sample  grid  has  188  cells  in  the  cross-shore  direction  and  236  cells  in  the 
alongshore  direction.  The  grid  spacing  is  25  m. 


Incident  Wave  Spectra 

Incident  two-dimensional  wave  spectra  are  specified  as  energy  density  as  a 
function  of  frequency  and  direction.  A  single  input  spectrum  is  applied  along  the 
entire  offshore  boundary  of  the  STWAVE  grid  (the  spectrum  is  set  to  zero  for 
any  land  points  along  the  boundary).  Thus,  it  is  good  practice  to  establish  the 
offshore  boundary  of  the  bathymetry  grid  along  a  depth  contour  where  the  wave 
spectrum  is  fairly  homogeneous  (no  large  shoals  or  canyons  offshore  of  the 
boundary).  The  first  line  of  the  input  spectra  file  describes  the  number  of 
frequency  and  direction  bins  in  the  spectra  (read  in  free  format): 

NF  =  Number  of  frequency  bins  in  the  spectra.  The  number  of  frequency 
bins  determines  how  finely  the  calculated  spectra  are  resolved.  A  large  number 
of  frequency  bins  increases  the  computation  time,  and  a  small  number  of  bins 
reduces  model  resolution.  Typically,  20-30  bins  are  used. 

NA  =  Number  of  direction  bins  in  the  spectra.  This  value  must  be  set  to  35, 
which  gives  5-deg  resolution  in  direction. 

The  next  lines  of  the  incident  wave  spectra  file  specify  the  frequencies  for 
model  spectra  (used  for  the  input  spectra,  internal  computations,  and  output 
spectra),  starting  from  the  lowest  frequency.  There  must  be  NF  frequencies 
specified.  The  frequencies  are  again  read  in  free  format  and  may  occupy  as 
many  lines  as  needed.  These  frequencies  should  span  the  frequency  range  where 
significant  wave  energy  is  contained  in  the  spectrum.  This  can  be  estimated  by 
inspecting  the  input  spectrum  or  estimating  the  peak  period  expected  using  the 
wave  growth  curves  in  the  Coastal  Engineering  Manual  (in  preparation).  A  rule 
of  thumb  is  that  the  spectral  peak  should  fall  at  about  the  lower  one-third  of  the 
frequency  range  (e.g.,  if  the  peak  frequency  is  0.1  Hz,  the  range  may  be  0.01  to 
0.3  Hz).  Wave  frequencies  higher  (or  periods  shorter)  than  the  highest  frequency 
bin  or  lower  than  the  lowest  frequency  bin  will  not  be  resolved  by  the  model. 
Typically,  frequency  increments  are  on  the  order  of  0.01  Hz,  but  the  increment 
need  not  be  constant.  West  Coast  applications  will  tend  to  require  finer 
resolution  focused  at  lower  frequencies  (because  of  long  wave  periods),  and  Gulf 
Coast  or  Great  Lakes  applications  will  tend  to  require  coarser  resolution  covering 
a  broader  range  of  frequencies  (because  of  shorter  wave  periods). 

Following  specification  of  the  frequency  bins  is  a  header  line  containing  a 
spectrum  identifier,  wind  information,  peak  frequency,  and  water  elevation 
correction: 
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IDD  =  Integer  event  identifier,  such  as  a  date  or  test  number.  The  integer  is 
limited  to  nine  digits. 

U  =  Wind  speed  in  meters/second.  This  wind  speed  is  considered  constant 
over  the  entire  grid.  Winds  that  are  blowing  offshore  (-x  direction)  are  neglected 
by  the  model,  because  the  model  only  treats  generation  in  a  half  plane.  For  wind 
directions  greater  than  60  deg  relative  the  x-axis,  STWAVE  underpredicts  wave 
generation  (because  of  the  half-plane  model  coverage).  To  model  large  ranges  of 
wind  direction,  multiple  STWAVE  grid  orientations  may  be  required. 

UDIR  =  wind  direction  relative  to  the  STWAVE  coordinate  system  in 
degrees.  Wind  direction  is  measured  counterclockwise  from  the  x-axis  (same 
convention  used  for  wave  direction,  see  Figure  3). 

FM  =  Peak  spectral  frequency  in  Hz.  If  a  spectrum  is  specified  with  no 
energy  (i.e.,  waves  will  be  locally  generated  in  the  modeling  domain),  FM  should 
be  set  to  the  highest  frequency  specified  in  the  model  parameter  file.  This 
parameter  is  no  longer  used  in  the  STWAVE  computations,  but  is  calculated 
from  the  input  spectrum.  A  value  for  FM  must  be  included  in  the  input  spectra 
file  to  maintain  backward  compatibility  with  STWAVE  version  2.0. 

DADD  =  Water  elevation  correction  (tide  or  storm  surge)  in  meters,  relative 
to  the  bathymetry  datum.  Users  should  take  special  care  to  ensure  a  consistent 
datum  is  used  to  specify  the  bathymetry  and  the  water  elevation  correction.  The 
water  elevation  correction  is  constant  over  the  entire  grid.  Positive  values  of 
DADD  increase  water  depth  and  negative  values  decrease  water  depths. 

The  header  line  is  read  in  free  format  and  is  followed  by  the  energy  densities 
in  the  units  of  meters  squared/hertz/radian.  The  spectrum  is  read  starting  with  the 
lowest  frequency  and  reading  all  the  wave  directions  (from  -85  deg  to  85  deg), 
then  reading  energy  densities  for  all  directions  for  the  next  lowest  frequency,  etc. 
The  total  number  of  energy  densities  read  equals  the  number  of  frequencies  times 
the  number  of  directions. 

Multiple  wave  model  runs  can  be  executed  by  inserting  multiple  spectra, 
with  identifying  header  lines,  in  the  input  file  (the  number  of  frequency  and 
direction  bins  and  the  frequencies  are  not  repeated  for  each  spectrum,  but  must 
remain  the  same).  The  model  will  loop  through  the  execution  until  it  runs  out  of 
input  spectral  information  (model  parameters  and  bathymetry  will  stay  constant). 
If  wave-current  interaction  is  specified  (ICUR=1),  the  number  of  current  fields  in 
the  current  input  files  must  equal  the  number  of  input  spectra.  The  default  name 
for  the  incident  spectra  file  is  spec .  in.  In  SMS,  the  file  is  specified  as 
pro  j  ect .  eng  (.eng  for  energy  density).  A  sample  incident  spectrum  file  is 
given  in  Appendix  C.  The  sample  spectrum  file  has  25  frequencies  and 
35  directions.  The  spectral  frequencies  are  from  0.01  to  0.25  Hz,  incrementing 
by  0.01  Hz.  The  event  identifier  is  200700  (12  March  1996  at  0200  GMT).  The 
wind  speed  is  15  m/sec,  wind  direction  is  10  deg,  and  the  tide/surge  elevation  is 
-1 .0  m  relative  to  mean  sea  level  (which  corresponds  to  the  mean  sea  level  datum 
used  to  develop  the  bathymetry  file). 
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The  incident  wave  spectrum  can  be  specified  using  a  measured  two- 
dimensional  spectrum,  by  generating  a  spectrum  using  a  parametric  shape,  or 
from  a  spectrum  calculated  by  a  global-  or  regional-scale  wave  model.  Rarely  do 
measured  spectra  have  the  directional  resolution  required  by  the  model;  thus, 
measured  spectra  require  interpolation  to  finer  resolution.  A  parametric  spectral 
shape  together  with  a  directional  spreading  function  can  also  be  applied  to 
specify  an  incident  spectrum,  knowing  the  wave  height,  period,  and  direction. 
SMS  includes  the  capability  to  generate  incident  spectra  using  a  TMA  one¬ 
dimensional  shallow-water  spectral  shape  (named  for  the  three  data  sets  used  to 
develop  the  spectrum:  TEXEL  storm,  MARSEN,  and  ARSLOE)  (Bouws  et  al. 
1985)  and  a  cos””a  directional  distribution  (see  Chapter  5).  To  generate  a  TMA 
spectrum,  the  following  parameters  must  be  specified:  peak  wave  period  (Tp), 
wave  height,  water  depth,  and  a  spectral  peakedness  parameter  (y).  The 
peakedness  parameter  controls  the  width  of  the  frequency  spectrum  (small 
numbers  give  broad  peaks  and  large  numbers  give  narrow  peaks).  The 
directional  distribution  of  the  spectrum  is  specified  with  a  mean  direction  and  a 
directional  spreading  coefficient  {mi).  The  energy  in  the  frequency  spectrum  is 
spread  proportional  to  cos nn(a-am),  where  a  is  direction  of  the  spectral 
component  and  am  is  the  mean  wave  direction.  Guidance  for  selecting  yand  nn 
is  given  in  Table  1  (Thompson  et  al.  1996).  Spectra  with  large  peak  periods  are 
generally  narrow  in  both  frequency  and  direction  (swell).  For  small  peak 
periods,  spectra  are  typically  broad  in  both  frequency  and  direction  (sea). 


Table  1 

Approximate  Spectral  Peakedness  and  Directional  Spreading 
Parameters 

Tp ,  Sec 

r 

nn 

<  10 

3.3 

4 

11 

4 

8 

12 

4 

10 

13 

5 

12 

14 

5 

16 

15 

6" 

18 

16 

6 

20 

17 

7 

22 

18 

7 

26 

19 

8 

28 

20 

8 

30 

Current  Fields 

A  current  field  input  file  is  required  only  if  ICUR  =  1  or  2  in  the  model 
parameter  file  (it  is  ignored  if  ICUR  =  0).  The  current  vector  is  specified  as 
x-  and  y-components  of  the  current  at  each  grid  cell,  u  and  v,  in  the  units 
meters/second.  The  x-  and  y-directions  correspond  to  the  x-  and  y-axes  of  the 
grid  in  the  local  STWAVE  coordinates  (Figure  3).  The  first  line  of  the  current 
input  file  is  NI,  NJ,  and  DXINC,  identical  to  the  parameters  specified  in  the 
bathymetry  file.  STWAVE  checks  that  these  values  are  consistent  with  the 
bathymetry  file,  and  if  they  are  not,  model  execution  is  terminated  with  the  error 
“current  field  does  not  match  depth  grid  size.”  Next,  the  header  for  the  current 
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field  is  given,  which  is  an  event  identifier.  This  header  is  not  used  by  the  model, 
but  is  helpful  to  ensure  that  the  current  and  spectral  input  files  coincide.  Similar 
to  the  bathymetry  file,  the  current  file  begins  with  the  u/v  pair  at  cell  (1,  NJ)  and 
reads  in  the  cross-shore  direction  (I  =  1  to  NI).  Then  the  read  is  repeated  for  J  — 
NJ-1  (I  =  1  to  NI),  and  progresses  to  J  =  1 .  The  current  components  are  read  in 
free  format.  The  default  name  for  the  current  field  file  is  current .  in.  In 
SMS,  the  file  is  specified  as  pro  j  ect .  cur.  A  sample  current  field  file  is  given 
in  Appendix  D.  In  the  sample  file,  the  first  line  repeats  the  grid  parameters  from 
the  bathymetry  file.  The  grid  is  defined  as  1 88  cells  in  the  cross-shore  direction 
and  236  cells  in  the  alongshore  direction  with  a  grid  spacing  of  25  m.  The  event 
identifier  for  the  current  field  is  200700,  which  matches  the  event  identifier  in  the 
incident  spectrum  file.  As  with  the  spectral  file,  multiple  wave  model  runs  can 
be  executed  by  inserting  multiple  current  fields  (with  their  event-identifier 
headers)  in  the  input  file  (as  well  as  the  same  number  of  input  spectra  in  the 
spec .  in  file)  for  ICUR  =  1 .  The  model  will  loop  through  the  execution  until  it 
runs  out  of  input  spectral  and  current  information.  For  ICUR  =  2,  one  current 
field  will  be  used  for  all  model  runs  (the  first  or  only  current  field  in  the 
current .  in  file).  This  is  a  new  option  in  STWAVE  version  3. 
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4  Output  File  Description 


STWAVE  has  up  to  five  output  files.  These  files  contain  wave  spectra  at  the 
selected  output  points;  wave  height,  period,  and  direction  at  the  selected  output 
points;  fields  of  wave  height,  period,  and  direction  over  the  entire  STWAVE 
modeling  domain;  fields  of  breaker  indices  indicating  active  regions  of  breaking 
over  the  entire  domain;  and  fields  of  radiation  stress  gradients  over  the  entire 
domain.  If  multiple  incident  wave  spectra  are  specified  in  the  input  file  (and 
current  fields,  if  wave-current  interaction  is  specified),  then  output  from  each  of 
these  model  runs  is  appended  to  the  output  files. 


Selected  Wave  Spectra 

Wave  spectra  at  selected  grid  cells,  specified  in  the  model  parameters  file 
(see  Chapter  3,  “Input  File  Description”),  are  saved  in  a  spectral  output  file. 
Similar  to  the  input  spectra  file,  the  first  line  of  the  output  file  describes  the 
number  of  frequencies  and  number  of  directions  (NF  and  NA),  and  the  following 
lines  give  the  frequencies.  Each  spectrum  in  the  file  is  proceeded  by  the 
following  header  information: 

IDD  =  Integer  event  identifier,  specified  in  the  incident  wave  spectrum  file. 

I  =  Integer  x  grid  cell  index  of  the  special  output  point,  specified  in  the  model 
parameter  file. 

J  =  Integer  y  grid  cell  index  of  the  special  output  point,  specified  in  the  model 
parameter  file. 

N  =  Integer  index  of  the  special  output  point  (values  are  1  to  NSELCT,  based 
on  the  order  in  which  the  points  are  specified  in  the  model  parameter  file). 

The  header  format  is  (il0,3i5).  Following  the  header  line,  the  output  spectrum  is 
written  in  the  same  order  and  units  in  which  the  incident  spectrum  was  read. 
Energy  densities  are  given  in  the  units  meters  squared/hertz/radian.  The 
spectrum  is  written  starting  with  the  lowest  frequency  and  writing  energy  density 
for  all  wave  directions  (from  -85  deg  to  85  deg),  followed  by  the  energy  densities 
for  all  directions  for  the  next  lowest  frequency,  etc.  The  format  of  the  energy 
densities  is  (17f8.3). 
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The  spectra  will  be  written  to  the  file  in  an  order  based  on  the  I  cell  index, 
starting  from  the  smallest  index  and  proceeding  to  the  largest  index.  This  is 
because  spectra  are  not  saved  over  the  entire  grid,  so  the  output  spectra  are 
written  to  the  file  as  the  grid  column  is  solved  and  then  are  overwritten  in 
memory.  Thus,  the  order  may  not  coincide  with  N=  1,  NSELCT.  The  default 
file  name  for  the  output  spectra  file  is  spec .  out.  In  SMS,  the  file  is  specified 
with  a  project  name  and  the  extension  .obs. 


Selected  Wave  Parameters 

The  wave  parameters  Hmo,  Tp,  and  mean  wave  direction,  am,  at  selected  grid 
cells,  specified  in  the  model  parameters  file  (see  Chapter  3,  “Input  File 
Description”),  are  saved  in  a  selected  parameter  output  file.  The  selected  grid 
cells  are  the  same  locations  where  spectra  are  saved.  For  each  output  point,  the 
following  information  is  written: 

IDD  =  Integer  event  identifier,  specified  in  the  incident  wave  spectrum  file. 

I  =  Integer  x  grid  cell  index  of  the  special  output  point,  specified  in  the  model 
parameter  file. 

J  =  Integer  y  grid  cell  index  of  the  special  output  point,  specified  in  the  model 
parameter  file. 

Hm0  =  Zero-moment  wave  height  at  cell  (I,  J),  in  meters. 

Tp  =  Peak  wave  period  at  cell  (I,  J),  in  seconds. 

am  =  Mean  wave  direction  at  cell  (I,  J),  in  degrees  relative  to  the  STWAVE 
grid  (see  Figure  3). 

The  format  of  each  line  is  (il0,2i3,f6.2,f6.1,f6.0).  Information  for  each  selected 
grid  cell  is  written  on  a  separate  line  in  the  file.  The  parameters  are  written  to  the 
file  in  an  order  based  on  the  I  cell  index,  starting  from  the  smallest  index  and 
proceeding  to  the  largest  index.  The  default  file  name  is  selhts.  This  file  is 
not  used  in  SMS  because  wave  height,  period,  and  direction  for  all  grid  cells  are 
displayed  graphically  on  the  screen  over  the  entire  domain  (values  at  any  cell  in 
the  domain  are  displayed  by  clicking  on  the  cell  with  the  mouse).  A  sample  of  a 
wave  parameter  file  is  given  in  Appendix  E. 


Wave  Parameter  Fields 

The  wave  parameters  Hmo,  Tp,  and  am  are  also  saved  for  all  grid  cells  into  a 
wave  field  file.  Written  first  in  the  file  are  the  grid  dimensions  (NI  and  NJ)  and 
grid  spacing  (DXINC)  repeated  from  the  bathymetry  input  file.  Next  in  the  file  is 
the  event  identifier  (IDD  from  the  incident  wave  spectrum  file)  in  the  format 
(ilO).  After  this  header  line,  all  the  wave  heights  are  written,  beginning  at  cell  (1, 
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NJ)  and  proceeding  in  the  cross-shore  direction  (1=1,  NI).  The  write  is  repeated 
for  each  grid  row  from  J  =  NJ-1  to  J  =  1 .  The  format  is  (1 6f6.2).  Following  the 
wave  heights,  the  wave  periods  are  written  in  the  same  order  with  the  format 
(16f5.1).  Then  the  wave  directions  are  written  (again  in  the  same  order)  with  the 
format  (16i5).  If  multiple  input  spectra  are  specified,  the  event  identifier  and  the 
wave  heights,  periods,  and  directions  for  each  condition  are  appended  to  the  end 
of  the  file. 

The  default  file  name  for  the  wave  field  file  is  wavf  Id.  In  SMS,  the  file  is 
specified  with  a  project  name  and  the  extension  .  wav. 


Breaker  Index  Fields 

The  regions  where  waves  are  actively  breaking  are  of  interest  for  sediment 
transport  and  navigation  studies.  Grid  cells  with  active  breaking  (see 
Equation  12)  are  specified  with  an  index  of  1  and  nonbreaking  cells  with  an 
index  of  0.  These  indices  for  the  entire  domain  are  written  to  the  file  break  (or 
pro  j  ect .  brk  in  SMS)  if  the  input  parameter  IBREAK  =  1  (see  Chapter  3). 
Written  first  in  the  file  are  the  grid  dimensions  (NI  and  NJ)  and  grid  spacing 
(DXINC)  repeated  from  the  bathymetry  input  file.  After  this  header  line,  all  the 
breaker  indices  are  written,  beginning  at  cell  (1 ,  NJ)  and  proceeding  in  the  cross¬ 
shore  direction  (1=1,  NI).  The  write  is  repeated  for  each  grid  row  from  J  =  NJ-1 
to  J  =  1 .  The  format  is  (1 6i5).  If  multiple  input  spectra  are  specified,  the  breaker 
indices  for  each  condition  are  appended  to  the  end  of  the  fi  le.  No  file  is  written 
when  the  input  parameter  IBREAK  =  0.  A  partial  listing  of  a  sample  wave 
breaking  file  is  given  in  Appendix  F. 


Radiation  Stress  Gradient  Fields 

The  radiation  stress  gradients  rjpw  and  t/pw  are  calculated  at  all  grid  cells 
when  IRS  =  1.  These  stress  gradient  fields  are  written  to  the  file  radstress. 
Written  first  in  the  file  are  the  grid  dimensions  (NI  and  NJ)  and  grid  spacing 
(DXINC)  repeated  from  the  bathymetry  input  file.  Next  in  the  file  is  the  event 
identifier  (IDD  from  the  incident  wave  spectrum  file)  in  the  format  (i  1 0) .  After 
this  header  line,  the  (r/pw,  r/p„)  pair  are  written  beginning  at  cell  (1,  NJ)  and 
proceeding  in  the  cross-shore  direction  (1=1,  NI).  The  write  is  repeated  for  each 
grid  row  from  J  =  NJ-1  to  J  =  1.  The  stress  gradients  are  written  in  free  format. 

If  multiple  input  spectra  are  specified,  the  event  identifier  and  the  radiation  stress 
gradients  for  each  condition  are  appended  to  the  end  of  the  file.  A  partial  listing 
of  a  sample  radiation  stress  gradient  file  is  given  in  Appendix  F. 
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5  SMS  User  Interface 


Information  contained  in  the  STWAVE  input  files  (Chapter  3)  and  output 
files  (Chapter  4)  is  conceptually  easy  to  understand,  but  practically,  these  files 
contain  large  amounts  of  data  that  can  be  difficult  to  generate  and  interpret.  To 
assist  STWAVE  users  in  generating  input  files  and  visualizing  output  files,  a  user 
interface  has  been  built  for  STWAVE  within  the  Surface-water  Modeling  System 
(Brigham  Young  University  Engineering  Computer  Graphics  Laboratory 
(ECGL)  1997).  The  SMS  interface  supports  grid  generation,  interpolation  of 
current  fields,  generation  of  input  spectra,  visualization  of  wave  heights,  periods, 
and  directions,  and  visualization  of  output  spectra.  U.S.  Army  Corps  of 
Engineers  users  are  supported  through  the  Coastal  and  Hydraulics  Laboratory 
('http://chl.wes.armv.mil/software/sms/).  and  non-Corps  users  can  get  information 
from  ECGL  (http://www.ecgl.bvu.edu'). 

STWAVE  is  a  self-contained  numerical  model  capable  of  executing  from  a 
DOS  command  prompt.  However,  SMS  has  a  comprehensive  environment  that 
provides  powerful  pre-  and  post-processing  tools  to  develop  and  visualize 
STWAVE  input  and  output  files  as  well  as  launch  an  external  execution  of 
STWAVE. 

This  chapter  provides  a  brief  tutorial  on  an  application  of  STWAVE  within 
SMS  version  7.1  Beta  (October  2000).  It  does  not  cover  all  the  SMS  features, 
but  describes  a  typical  application  of  the  model  starting  with  grid  development 
from  bathymetry,  generation  of  multiple  input  spectra,  monitoring  station 
identification,  interpolation  of  multiple  current  fields  from  an  ADCIRC  finite- 
element  bathymetry  grid,  model  control,  model  execution,  and  visualization  of 
the  model  results. 

Commands  in  SMS  are  accessed  through  pulldown  menus  located  on  the 
menu  bar  (Figure  5)  along  the  top  of  the  Graphical  User  Interface  (GUI) 

(Figure  6).  The  edit  window  (Figure  5)  lies  below  the  menu  bar.  It  includes 
controls  editing  coordinates  (type  in  a  new  value  and  press  the  enter  or  tab  key). 
Powerful  modules,  tools,  and  macros  are  accessible  from  icons  in  the  tool  palette 
(Figure  6)  on  the  left  side  of  the  GUI  below  the  edit  window. 

The  SMS  GUI  supports  six  distinct  modules,  one  for  each  of  the  following 
datatypes:  mesh  SL  grid  ffl,  Cartesian  IE.  scatterpoint  El,  map  ffl,  and  river 
modules  HI  The  six  module  icons  are  located  at  the  top  of  the  tool  palette.  As 
you  switch  from  one  module  to  another,  the  tool  palette  and  menu  bar  changes. 
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File  Edit  Display  Data  5TWAVE  Window  Help 
Coverage:  :  default  coverage  ^  Solution:  :  Generic  Datasets 


X  Location:  50.30 


V  Location:  508.14 


▼1  Time  Step:  i  No  Til 


Z  Location:  ? 


Figure  5.  SMS  GUI  menu  bar  and  edit  window 


File  Edit  Display  Data  STWAVE  Window  Help 
Coverage:  default  coverage  y  |  Solution:  Generic  Datasets 


X  Location:  2.96 


Y  Location:  411.98 


Time  Step:  No  Time 


Z  Location: 


Figure  6.  SMS  GUI 


This  allows  the  user  to  focus  only  on  the  tools  and  commands  related  to  the  data 
you  are  currently  processing. 

To  promote  clarity  throughout  this  chapter,  the  SMS  GUI  menu  selections 
are  expressed  as  Menu  |  Selection  and  modules  are  expressed  as  DataType 
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Module.  SMS  saves  all  its  working  data  and  settings  in  project  files.  These 
SMS  created  project  files  are  given  the  project’s  name  followed  by  an  extension. 
In  this  chapter,  the  primary  SMS  project  file  is  expressed  as  pro  j  ect .  spr 
while  other  SMS  project  files  are  designated  in  the  text  as  project .  ext. 

Build  STWAVE  Cartesian  Grid 

To  begin  building  a  grid  in  SMS,  the  input  bathymetry  (xyz)  file  should  be  in 
column  format,  with  comma-,  tab-,  or  space-delimited  fields.  The  bathymetry 
file  may  include  a  file  header  of  any  number  of  lines.  These  lines  may  be  used 
for  recording  the  coordinate  system,  vertical  datum,  source  of  data,  point  of 
contact,  and  other  pertinent  information.  SMS  can  perform  coordinate  system 
conversions,  add  information  from  nautical  charts  by  digitizing  sounding  depths, 
and  import  map  images. 

The  following  series  of  SMS  instructions  illustrates  the  development  of  a  2-D 
grid  file  (pro  j  ect .  dep)  that  is  used  as  the  input  bathymetry  file  for  STWAVE 
(Chapter  3,  “Bathymetry”): 

Start  SMS. 

a.  Open  and  read  bathymetry  file,  YAQbath.xyz,  and  make  necessary 
coordinate  system  conversions: 

1 .  Scatter  Module  -  E 

(a)  File  |  Open->YAQbath.xyz. 

(1)  File  Import  Wizard  dialog  automatically  displays.  Apply 
defaults  with  these  exceptions: 

i.  Scatter  Options:  Tolerance  =  .1 

ii.  Convert  bathymetry  units  from  feet  to  meters. 

♦  Current  Horizontal  System: 

■  State  Plane  NAD27,  Oregon  North,  meters. 

■  NGVD27,  meters. 

♦  Convert  incoming  data  from: 

■  State  Plane  NAD27,  Oregon  North,  feet. 

■  NGVD27,  feet. 
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Figure  7.  Bathymetry  scatter  set  with  excess  river  points  featured  in  black  area 


b.  Remove  excess  river  points  from  bathymetry  scatter  set.  River  points  to 
be  removed  are  outlined  in  black  in  Figure  7. 

1 .  Scatter  Module  -  E3 

(a)  Edit  |  Select  with  poly... 

(1)  Click  mouse  at  points  around  river  section  to  be  removed. 

(2)  Double-click  to  close  selection  area. 

(3)  Edit  |  Delete  to  remove  points. 

(4)  Repeat  the  previous  steps  as  needed. 

c.  STWAVE  requires  positive  Z-coordinate  values.  YAQbath.xyz 
bathymetry  file  has  negative  Z-coordinates  that  SMS  can  convert  to 
positive  Z-coordinates  by: 

1.  Scatter  Module -0 
(a)  Data  |  Transform- 

(1)  Depth<->Elevations. 

d.  The  land/ocean  boundary  can  be  improved  for  this  data  set  by  importing  a 
scanned  image  (18581a .  tif)  of  the  Yaquina  area  from  NOS  chart 
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1 8581  to  display  in  SMS  with  the  bathymetry  points  (aerial  photographs 
may  also  be  used).  The  image  should  be  in  TIFF  format. 

1 .  Scatter  Module  -  L±. 

(a)  File  |  Open ->  18581a.tif.  The  Register  Image  dialog 
automatically  appears.  See  the  SMS  “Register  Image”  Help  topic 
for  further  explanation  of  registering  an  image. 

(b)  To  register  the  image  for  use  as  a  background  display,  you  need 
to  know  the  coordinates  of  three  locations  (street,  railroad,  or 
latitude/longitude  point  on  a  map,  etc.)  on  the  image.  Align  each 
+  symbol  with  its  known  reference  point  on  the  image.  This 
alignment  defines  the  three  U,V  coordinates  (local  image 
coordinates)  in  the  edit  fields  below  the  image.  Then,  enter  the 
X,Y  values,  which  are  the  state  plane  coordinates  in  meters 
corresponding  to  these  reference  points,  in  the  edit  fields  below 
the  image. 

(1)  To  enlarge  an  area  on  the  image,  click  the  Zoom  Tool  E 
and  drag  the  mouse  cursor  over  the  +  symbol  to  form  a 
rectangle  (click  the  Frame  Macro  Tool  |Vl  to  zoom  out). 
For  this  example,  the  user  can  directly  input  the  U,  V,  X, 
and  Y  to  properly  register  the  image. 

i.  Point  #1 

1.  U:  660 

2.  V:  2134 

3.  X:  329335.841 

4.  Y:  113546.3735 

ii.  Point  #2 

1.  U:  622 

2.  V:  578 

3.  X:  325372.4217 

4.  Y:  113721.0504 

iii.  Point  #3 

1.  U:  2807 

2.  V:  525 

3.  X:  325126.0839 

4.  Y:  108170.6839 

(c)  Click  OK. 
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e.  Land  boundaries  are  not  defined  in  the  scatter  set.  The  following 

instructions  set  the  shoreline  to  0.0  m  and  the  land  to  -2.0  m  by  digitizing 
those  features  from  the  NOS  chart  image: 

1.  Map  Module -S 

(a)  FeatureObjects  |  Coverage 

(1)  CoverageType ->  STWAVE. 

(b)  ToolPallette->Create  Feature  Arc  Tool  -  E 

(c)  Click  mouse  cursor  along  shoreline.  Continue  clicking  the  cursor 
around  the  land  area  until  the  arc  is  completed. 

(d)  FeatureObjects  |  BuildPolygon.  In  response  to  the  dialog 
asking  whether  to  Use  all  arcs?,  click  OK. 

(e)  Select  the  Feature  Polygon  Tool  -  IM  Then,  double-click 
inside  the  land  polygon.  Select  Land  and  OK. 

f  To  edit  or  delete  a  node: 

1.  Map  Module -ff 

(a)  FeatureObjects  |  Coverage  ->  CoverageType  -> 
STWAVE. 

(b)  ToolPallette->Select  Feature  Arc  Tool  -  E 

(c)  Click  on  the  node  -> 

(1)  hold  down  left  mouse  button  ->  drag  node  to  reposition. 

(2)  Select  Edit  |  Delete  to  delete  the  node. 

g.  Mapping  to  the  Grid  defines  the  grid  origin,  orientation,  and  cell  size: 

1.  Map  Module  -K 

(a)  Feature  Objects  ->Grid  Frame->New  Grid: 

(1)  Origin:  x=321470  m,  y=108765  m. 

(2)  Orientation:  350  deg.  The  orientation  is  the  direction  of 
the  x-axis,  measured  counterclockwise  from  East.  This 
information  is  used  in  SMS  to  visualize  the  grid  in  the 
correct  orientation. 
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(3)  Grid  size:  4,700  m  by  5,890  m. 

(4)  OK 


30 


(b)  Feature  Objects  ->  Map  ->  2D  Grid 

(1)  Cell  size  =  25  m 

(2)  Depth  =  Interpolated  (use  defaults). 

(3)  OK 

2.  Alternatively,  the  grid  position  may  be  defined  graphically.  Rotate 
the  grid  by  putting  the  cursor  on  the  open  circle  in  the  lower  right 
comer  of  file  grid  and  dragging  to  a  new  orientation.  Translate  the 
grid  by  grabbing  the  purple  frame  or  resize  by  grabbing  the  black 
circles  on  the  grid  comers  or  sides.  The  grid  origin  is  the  comer  with 
the  I,J  axis.  Refer  to  the  SMS  Grid  Frame  help  topic  for  more 
information. 


h.  Map  jetty  locations.  Jetty  elevations  are  not  included  in  the  xyz 

bathymetry,  so  the  jetty  elevations  must  be  added  to  the  grid.  A  scanned 
image  (yaq_entrance.tif)  of  the  jetty  location  in  Tiff  format  is 
imported  to  visualize  the  jetty  location.  This  image  (yaq_entrance.tif) 
is  more  current  than  the  NOS  image  (1 8581  a.tif): 


1. 


Map  Module  - 
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(a)  File  |  Open  ->  yaq_entrance.tif. 

This  image  will  replace  any  previous  image  ->  OK. 

(b)  Register  the  image  for  use  as  the  grid’s  background  display  by 
identifying  three  points  on  the  image  (U,V)  which  correspond  to 
locations  having  known  real  world  (X,Y)  coordinates. 


(1)  Coordinates  of  the  three  registration  points,  shown  by  a  + 
symbol,  are  listed  in  edit  fields  below  the  image.  Drag  the 
large  +  symbol  over  the  small  +  symbol  having  a 
corresponding  number.  Click  the  Zoom  Tool  cSJ  and  drag 
the  mouse  cursor  over  the  +  symbol  to  form  a  rectangle. 
This  enlarges  the  area  inside  the  rectangle.  Ensure  that  the 
large  &  small  +  symbols  are  aligned.  If  not,  click  on  the 
Select  Point  Tool  IS  and  adjust  the  small  +  symbol  so 
that  it  is  centered  over  the  large.  Click  the  Frame  Macro 
Tool  S  to  zoom  out. 


(2)  Enter  State  Plane  coordinates  (in  meters)  for  the  three 
registration  points. 
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i.  Point  #1: 

1.  X:  323798.2 

2.  Y:  109018.4 

ii.  Point  #2: 

1.  X:  324897.9 

2.  Y:  113080.5 

iii.  Point  #3: 

1.  X:  326104.3 

2.  Y:  110565.9 

iv.  OK 


i.  Turn  off  bathymetry  scatter  points: 
1.  Scatter  Module -1 


F* 

4- 


(a)  Display  |  Display  Option. 

(b)  Remove  check  from  Points  by  clicking  the  checked  box. 

j.  Change  elevations  along  the  jetty  to  -4.8  m,  which  is  the  jetty  crest 
elevation: 


1.  Cartesian  Grid  Module  -  H 


(a)  Display  |  Display  Options. 

(1)  Check  Cells. 

(2)  Press  Contours  Options  button. 

i.  Contour  Method:  Normal  Linear 

ii.  OK 

(3)  OK. 

(b)  Select  the  Zoom  Tool  0.  Drag  the  mouse  cursor  over  the  jetty 
to  zoom  in  on  the  jetty. 

(c)  Select  the  Cell  Tool  H.  Hold  down  shift  key  and  click  the 
mouse  on  cells  to  be  modified.  Select  the  jetty  grid  cell  and  at 
least  one  grid  cell  on  either  side  of  the  jetty  center  line  to  total 

3  grid  cells.  Grid  cells  selected  will  be  color  enhanced.  From  the 
Edit  Window,  change  the  Z  coordinate  to  -4.8  and  press  the 
Enter  key.  All  selected  cells  are  changed  to  -4.8-m  elevation. 

(d)  Zoom  out  by  pressing  Frame  the  Data  Set  [v|  in  the  Tool 
Palette. 


Chapter  5  SMS  User  Interface 


31 


(e)  Cartesian  grid  completed. 
k.  Turn  off  grid  frame,  nodes  &  arcs: 


1. 


Map  Module  - 


m 


(a)  Display  |  Display  Options  ->  Remove  check  from  Nodes, 
Arcs,  and  Grid  Frame. 


/.  File  |  Save  Project...  Saves  all  the  SMS  working  files,  updates 
pro  j  ect .  spr,  and  creates  the  STWAVE  2-D  input  bathymetry  grid 
file  pro  j  ect .  dep .  Grid  building  may  be  temporarily  terminated  at 
any  time  in  the  process  previously  outlined  and  resumed  at  a  later  time  if 
the  project  was  saved. 

m.  To  resume  the  project: 

1.  File  |  Open ->  project.spr 


Generate  Input  Spectra  for  STWAVE 

SMS  can  generate  one  or  more  2-D  energy  spectra  and  save  the  spectra  in  a 
file,  projecteng.  Parameters,  which  are  entered  in  the  Spectral  Energy  dialog  for 
generating  spectra,  are  saved  by  SMS  in  pro  j  ect .  txt.  This  file, 
pro  j  ect .  txt,  is  not  required  by  SMS  or  STWAVE.  If  the  file  is  available, 
SMS  will  display  the  contents  in  the  Spectral  Energy  dialog  providing  you  with 
parameters  used  for  generating  the  2-D  energy  spectra. 

The  number  of  input  wave  spectral  events  determines  the  number  of 
STWAVE  model  simulations  performed.  This  example  features  two  events 
(spectra)  having  identifiers  175500  and  200700.  When  running  multiple 
simulations,  the  number  of  events  and  the  event  identifiers  should  be  the  same  for 
spectra  and  current  fields  (Chapter  3,  “Incident  Wave  Spectra”).  For  this 
example,  5-  and  10-m,  15-sec  waves  from  the  west  will  be  generated  with  ebb 
and  flood  current  fields. 

Begin  SMS. 

a.  Cartesian  Grid  Module  -  US 

1.  Stwave  |  Spectral  Energy ->  Solution=Generic  Data  Sets 

(Default) 

2.  Select,  Generate  Spectra. 

3.  Enter  input  spectral  parameters.  Guidance  on  selecting  these 
parameters  is  given  in  Chapter  3. 
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(a)  F  num  =  25  (25  frequencies) 

(b)  F  min  =  0.01  (Hz) 

(c)  F  delta  =  0.01  (Hz) 

(d)  Water  depth  =  47  (Depth  at  offshore  boundary  in  meters) 

(e)  Wave  dir  =  Global  (Use  standard  wave  direction  convention  or 
select  Local  to  use  STWAVE  local  convention  as  defined  in 
Figure  3.) 

(f)  Check  option  to  Replace  Old  Spectra. 

(g)  Number  of  Rows:  2  (number  of  events  to  generate) 

(h)  Index  =  175500,  H  =  5  m,  T  =  15  sec,  Angle  =  250  deg, 
WindSpeed  =15  m/sec,  WindDir  =  270  deg,  Tide  =  -1.0  m. 

(i)  Index=200700,  H=10  m,  T=15  sec,  Angle  =  250  deg, 
WindSpeed=15  m/sec,  WindDir=270  deg,  Tide  =  -1.0  m. 

(j)  Press  Generate. 

4.  If  spectra  are  not  visualized  on  the  Spectral  Energy  Dialog,  move 
the  cursor  around  the  window  and  look  at  the  top  of  the  dialog  for 
numbers  to  display  after  each  Cursor  parameter:  Freq,  Angle, 
and  Energy.  If  Cursor  parameter  values  do  not  appear,  this 
indicates  that  spectra  have  not  been  generated.  Repeat  steps  provided 
previously  to  generate  spectra.  If  Cursor  parameter  values  do 
appear,  click  the  Contour  Options  El  in  the  upper  left  comer  of 
the  Spectral  Energy  Dialog. 

(a)  Change  the  option,  Contour  between  Specified  Range,  to 

(1)  Maximum  Value  =  1 .0  (energy  density) 

(2)  Press  OK 

5.  To  plot  frequency  and  directional  distributions  of  spectra: 

(a)  Select  Freq.  Plot  button. 

(b)  Select  Dir.  Plot  button. 


6.  When  multiple  input  spectra  are  available  for  visualization,  choose 
one  by  selecting  a  time  or  integer  from  the  Time  Step  option  at  top 
right  of  the  Spectral  Energy  Dialog  or  click  the  Data  Browser 


and  select  a  Time  Step. 
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7.  Press  OK  to  close  the  Spectral  Energy  Dialog. 

b.  File  |  Save  Project...  Saves  all  the  SMS  working  files,  updates 
pro  j  ect .  spr,  and  creates  the  STWAVE  2-D  energy  spectra  file 
pro  j  ect .  eng.  If  the  project  was  saved,  development  of  input  spectra 
may  be  temporarily  discontinued  at  any  time  in  the  process  outlined 
above  and  resumed  at  a  later  time. 

c.  To  resume  a  project: 

1.  File  |  Open ->  project.spr 

Identify  Monitoring  Stations 

A  convenient  point  and  click  capability  is  available  in  SMS  for  the  user  to 
select  monitoring  stations.  Then,  monitoring  station  locations  are  saved  in  the 
STWAVE  input  file  pro  j  ect .  std  as  I-  and  J-grid-cell  indices  (Chapter  3, 
“Model  Parameters”).  After  running  the  model,  wave  spectra  for  each 
monitoring  station  are  written  to  the  STWAVE  output  file  pro  j  ect .  obs 
(Chapter  4,  “Selected  Wave  Spectra”).  Also,  model  results  for  wave  parameters 
Hmo,  Tp,  and  mean  wave  direction,  am,  are  written  to  the  STWAVE  output  file 
selhts  for  each  monitoring  station. 

Start  SMS 

a.  Cartesian  Grid  Module  -  US 

1 .  Display  monitoring  station  (cell)  on  the  grid  by  denoting  with  a  red 
box  inside  the  cell.  The  red  box  symbol  can  be  changed  in  the 
STWAVE  Display  Options. 

(a)  Display  |  Display  Options->check  Cells. 

(b)  Press  Options... 

(c)  Only  check.  Observation  Symbol  ->  OK 

(d)  OK 

2.  Repeat  the  following  steps  for  each  cell  you  wish  to  identify  as  a 
monitoring  station.  In  this  example  three  monitoring  stations  are 
identified  at  cell  locations:  (1 10,157),  (131,1 17),  (146,86) 

(a)  Click  Select  Cell  S  in  the  Tool  Palette. 

(1 )  OPTION:  Click  your  mouse  cursor  on  the  screen.  The  I,J 
location  of  the  cell  is  displayed  in  the  bottom  right  of  the 
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screen  in  the  edit  window.  Hold  down  the  shift  key  to 
select  more  than  one  cell. 

(b)  From  the  Menu  Bar,  select  DATA  |  Find  Cell 

(1)  You  may  enter  the  I,J  location  or  the  X,Y  coordinates  to 
identify  the  cell. 

(2)  Press  OK. 

(c)  STWAVE  |  Assign  Cell  Atts...  Check,  Monitoring  Station. 
Then,  press  OK. 

File  i  Save  Project...  Saves  all  the  SMS  working  files,  updates 
pro  j  ect .  spr,  and  creates  the  monitoring  station  file  proj  ect .  std. 


Wave-Current  Interaction 

Wave-current  interaction  requires  current  fields  imported  from  a 
hydrodynamic  model  (e.g.,  ADCIRC,  RMA2,  or  some  other  source)  on  the 
STWAVE  grid.  The  extent  of  the  current  model  domain  must  be  the  same  or 
larger  than  the  wave  model  domain.  Many  STWAVE  simulations  can  be  made 
using  one  current  field  (ICUR=2,  Chapter  3)  or  current  fields  may  be  specified 
for  each  input  wave  spectrum  (ICUR=1,  Chapter  3).  If  the  second  option  is 
chosen,  the  number  of  current  fields  must  equal  the  number  of  input  spectra.  For 
this  example,  there  are  two  input  current  fields  having  event  identifiers,  175500 
and  200700.  Because  wave-current  interaction  is  also  specified  for  this  example, 
the  number  of  input  spectra  events  equals  the  number  of  current  field  events,  and 
the  event  identifiers  for  input  spectra  match  the  event  identifiers  for  current  fields 
(Chapter  3,  “Current  Fields”). 

Two  ADCIRC  model  files  are  used  to  build  the  SMS  current  grid.  The  first 
file,  grid.  14,  is  an  ADCIRC  input  file  containing  the  finite-element  grid, 
bathymetry,  and  boundaries.  The  second  file,  Current_base.64,  is  an  ADCIRC 
output  file  containing  velocities  (u  and  v)  at  each  grid  node. 

The  following  SMS  instructions  create  a  2-D  current  grid  file 
(proj  ect .  cur)  that  is  used  as  the  input  current  field  file  for  STWAVE.  Note: 
create  a  new  directory  for  processing  current  field  data,  and  place  a  copy  of  files 
grid.  14  and  Current_base.64  in  the  new  directory.  This  precautionary  measure  is 
to  avoid  possibly  overwriting  the  STWAVE  Cartesian  grid  or  project  working 
files  outlined  in  “Build  STWAVE  Cartesian  Grid.” 

Start  SMS. 

a.  Mesh  Module  -  I® 

1.  Display  |  DisplayOption 
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(a)  Check:  elements,  contours,  and  mesh  boundary. 

2.  File  |  Open  ->  grid.14  as  a  file  of  type:  (*.grd,  *.14). 

3.  Data  |  DataBrowser  |  Import  ->  currentbase. 64  as  ADCIRC 
unit  64:  (*.64,  *.sol). 

(a)  Create  a  new  solution  set?  OK. 

(b)  Press,  Vector  ->  Scalars: 


(1)  Check,  All  Time  Steps 

(2)  Select  Vx  &  Vy. 

(3)  Perform  the  conversion. 

(4)  Done. 

(c)  Done. 

4.  Data  |  MeshToScatterPoint->scatterl,OK 
b.  Scatter  Module  - 1 


♦ 


1.  Display  |  DisplayOptions->Points  off. 

c.  Reduce  image.  The  area  to  retain  is  on  the  right  side  of  the  mesh 
approximately  halfway  down  where  you  will  notice  a  sinuous  feature. 
This  feature  is  the  river  referenced  previously  in  this  chapter  when 
discussing  how  to  “Build  STWAVE  Cartesian  Grid.”  Several  iterations 
of  the  following  commands  are  needed  to  remove  excess  points. 

1.  Mesh  Module-® 


(a)  Edit  |  Select  with  Poly... 

(1)  Click  mouse  around  the  grid  at  points  to  be  removed. 

(2)  Double-click  to  close  selection  area. 

(3)  Edit->Delete. 


(4)  Tool  Palette  |  Frame  Tool  t2l-  to  zoom  out. 

(5)  Repeat  1-4  as  needed. 

d.  Change  the  coordinate  system  of  the  input  current  field  to  the  coordinate 
system  of  the  wave  model. 
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1.  Scatter  Module  -l? 


(a)  Edit  |  CoordinateConversions: 

(1)  Select,  Current  Options  (for  incoming  data): 

i.  Horizontal=Geographic  NAD27 

ii.  Vertical=NGVD  29,  meters 

(2)  Select,  Convert  to  (same  coordinates  as  STWAVE  grid): 

i.  Horizontal=State  Plane  NAD  27,  Oregon  North,  meters 

ii.  Vertical=NGVD29,  meters 

(3)  Select,  Convert 

i.  A  dialog  may  appear  stating:  Unable  to  find  NAD27sp. 
Do  you  want  to  search  for  the  file? 

ii.  Select,  OK. 

iii.  NAD27sp  is  usually  in  directory 
C:\ProgramFiles\SMS70\latlon 

e.  Mapping  to  the  grid  defines  the  grid  origin,  orientation,  and  cell  size  to 
match  the  STWAVE  grid  previously  generated: 

1.  Map  Module -SI 

(a)  FeatureObjects  |  Coverages->STWAVE,  OK 

(b)  FeatureObjects  |  Grid  Frame 

(1)  X  =  321470  (m) 

(2)  Y=  108765  (m) 

(3)  Angle  =  350  (deg) 

(4)  Grid  size:  4700  (m)  by  5890  (m). 

(5)  Cols/Rows:  188,236 

(c)  FeatureObjects  |  Map->2-D  Grid 
(1)  Cell  size  =  25  (m) 
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(2)  Select,  Depth  Interpolated 

(3)  Press,  elevation. 

i.  Interpolated  Method:  Linear 

ii.  Scatter  Set:  scatter  1  (Same  name  assigned  in  prior 
mesh  to  scatter  set  conversion.) 

iii.  Function:  Magnitude  (64) 

iv.  Time-Steps:  Single 

v.  Time-Step  Interpolation:  Match  Time-Steps 

vi.  OK 

vii. NOTE:  This  example  is  using  high-resolution 
bathymetry  data  to  generate  the  STWAVE  grid.  When 
that  is  not  available,  this  step  can  also  generate  an 
STWAVE  depth  grid  based  on  the  ADCIRC 
bathymetry. 

(4)  Check,  Current. 

(5)  Press  Interpolated  to  interpolate  the  current  fields. 

(6)  Press,  Velocity  (64). 

i.  Interpolated  Method:  Linear. 

ii.  Scatter  Set:  Same  name  assigned  in  Depth 
Interpolated. 

iii.  Time-steps:  Multiple. 

iv.  Select,  Velocity64 

v.  Select,  X/Y 

vi.  Time-Step  Interpolation:  Match  Time-Step 

vii.  Beg=l  75500,  End=200700 

viii.  OK 

(7)  OK 

f  Turn  off  Mesh  Points: 
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1.  Mesh  Module  -  ffl 

(a)  Display  |  DisplayOptions 

(1)  All  off 

(2)  OK 

g.  Save  all  the  SMS  working  files,  update  proj  ect .  spr,  and  create  the 
STWAVE  2-D  input  current  field  file  proj  ec  t .  cur.  If  the  project  is 
saved,  current  field  file  building  for  input  to  STWAVE  may  be 
temporarily  discontinued  at  any  time  in  the  process  outlined  above  and 
resumed  at  a  later  time. 

1.  File  |  Save  Project... 

h.  To  resume  a  project: 

1.  File  |  Open ->  project.spr 

i.  Ensure  that  proj  ect .  cur  is  copied  to  the  STWAVE  project  directory, 
and  that  proj  ect .  cur  reflects  the  STWAVE  directory’s  project  name. 


Model  Control 

The  Model  Control  dialog  lets  the  user  select  appropriate  settings  for  the 
application  of  the  STWAVE  model.  These  settings  are  saved  by  SMS  in  the 
STWAVE  input  file,  project.std  (Chapter  3,  “Model  Parameters”). 

Start  SMS 

a.  Cartesian  Grid  Module  -  S 
1.  STWAVE  |  Model  Control... 

(a)  Check  grid  definition.  The  cell  size  and  number  of  columns/rows 
are  read  from  project.dep.  Do  not  change  these  values.  The 
origin  and  azimuth  values  (project.sim)  do  not  affect 
computations,  but  are  used  in  SMS  visualization. 

(1)  Cell  size  (m):  25. 

(2)  Number  of  columns:  188 

(3)  Number  of  rows:  236 

(4)  X  origin  (m):  321470. 
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(5)  Y  origin  (m):  108765. 

(6)  Azimuth  (deg):  350. 

(b)  Index:  Time-Step 

(c)  Source  Terms:  Source  Term  &  Propagation 

(d)  Wave  Current  Interaction:  Currents 

(e)  Breaking:  Write  Indices 

(I)  Radiation  Stresses:  Calculate 
(g)  OK 

b.  File  |  Save  Project. . .  Saves  all  the  SMS  working  files,  updates  project.spr, 
and  creates  the  STWAVE  input  parameter  file  project.std  (Chapter  3, 
“Model  Parameters”). 


Run  STWAVE 

When  all  the  necessary  input  files  are  created  (bathymetry,  input  spectra, 
monitoring  stations,  current  (optional),  and  model  control),  the  model  is  ready  to 
simulate  conditions  specified  by  the  input.  STWAVE  results  are  saved  in 
project. obs,  project.wav,  selhts,  project.brk,  and 
rads  tress  (Chapter  4) 

Start  SMS 

a.  Cartesian  Grid  Module  -B 
1.  File  |  Open ->  project.sim 
(a)  STWAVE  |  Run  model. 

( 1 )  A  dialog  appears  identifying  the  path  to  the  STWAVE 
executable. 

i.  Click  OK. 

(2)  Then,  a  DOS  window  is  launched  with  model  output 
directed  to  the  screen.  When  the  STWAVE  model  run 
completes,  the  message.  Press  return  to  continue,  is 
written  to  the  DOS  window.  The  DOS  window  disappears 
after  you  press  the  return  key. 
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STWAVE  may  be  ran  outside  of  SMS  from  the  command  prompt  by  entering  the 
command  stwave3  proj  ect .  sim.  Here,  stwave3  represents  the  model’s 
executable  filename  and  proj  ect .  sim  identifies  the  name  of  the  simulation 
file  saved  by  SMS.  If  the  user  wishes  to  build  a  proj  ect .  s  im  file  for  input 
into  SMS,  the  file  should  be  in  column  format,  with  comma-,  tab-,  or  space- 
delimited  fields. 

a.  The  following  example  proj  ect .  sim  file  specifies  the  STWAVE  grid 
origin,  orientation  azimuth,  and  the  STWAVE  input  and  output  files: 

STWAVE  xOrg  yOrg  azimuth 

DEP  pro j  ec  t . dep 

OPTS  project. std 

CURR  proj ect. cur 

SPEC  proj ect. eng 

WAVE  proj  ect . wav 

OBSE  proj ect. obs 

BREAK  proj ect. brk 

SPGEN  project.txt 


Load  Simulation 

Visualizing  STWAVE  model  results  is  possible  for  simulations  launched 
from  within  SMS  or  from  the  command  prompt  as  described  previously  in  “Run 
STWAVE.”  Because  STWAVE  is  a  self-contained  model,  its  model  results  are 
not  saved  within  SMS.  So,  when  a  simulation  finishes,  you  will  need  to  load  the 
simulation  file,  proj  ect .  sim,  into  SMS. 

The  directory  where  project.sim  resides  should  also  contain  the  files 
referenced  by  proj  ect .  sim  with  some  exceptions.  Three  optional 
proj  ect .  sim  files,  proj  ect .  cur,  proj  ect  .brk,  and  proj  ect .  txt 
may  not  be  in  the  directory.  If  the  optional  files  are  not  in  the  directory,  a 
warning  dialog  will  display  a  message  to  that  effect.  Click  OK  to  continue.  This 
helpful  warning  is  only  to  inform  the  user  of  file  status.  For  a  list  of  the  files 
referenced  by  project.sim,  see  the  previous  section,  “Run  STWAVE.” 

Start  SMS. 

a.  File  |  Open->  project.sim. 
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Visualize  Simulation  Output 

As  part  of  the  SMS  comprehensive  environment,  powerful  post-processing 
tools  are  provided  to  visualize  wave  information  estimated  by  the  model. 
Estimated  wave  information  is  identified  in  SMS  as  a  data  set.  Data  sets  consist 
of  a  scalar  or  vector  value  for  each  cell  in  the  grid.  Data  sets  are  considered 
steady  state  if  they  have  no  time  steps  and  are  considered  dynamic  if  they 
represent  different  points  in  time.  Refer  to  the  SMS  Help  topic  “Data  Sets”  for 
more  information. 

Use  the  Edit  Window  (Figure  8)  or  the  Data  Browser  to  select  Scalar  or 
Vector  data  set  information.  If  data  sets  are  representing  different  points  in 
time,  the  time-steps  for  the  data  sets  are  listed  in  the  Time-Step  list  box.  The 
identified  data  set  is  called  the  Active  Data  set. 


|SMS7. 

■■■SI:/ 

Jnlxi 

Fife  §* 

Display  Data  STWWE 

into  Help 

:c:- •••  -c 

Coverage; 

| default  coverage 

Sobta:  jbeakfSMAVE)  j 

2 

Soak  |Hei§li|  *|  4 

fete  [wave  | 

Y;:«c^r  . 

Zcoortiate  [26.822454 

S:  15?  VV  10.56 

%  0369 

X  Location: 

328823.73 

Y  Location:  114906.63 
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$;,?  ',V;  %  ? 

Vy  ? 

Figure  8.  SMS  edit  window 

Wave  height  estimates  are  displayed  as  the  default  SMS  scalar  data  set.  To 
view  wave  period  and  wave  direction  scalar  data  sets,  locate  the  Edit  Window  | 
Scalar  option  and  change  the  selection  to  period  or  direction.  You  may  also 
select  a  specific  time-step  from  the  EditWindow  |  TimeStep  list  box. 

Bathymetry  contours  may  be  superimposed  over  the  modeled  wave  field  to 
visualize  a  depth  frame  of  reference,  as  described  below. 


Visualize  and  store  bathymetry 

To  apply  bathymetry  contours  over  estimated  wave  information,  store 
bathymetry  contours  in  a  DXF  file.  Then,  load  the  DXF  file  to  have  the 
bathymetry  superimposed  on  contours  of  the  wave  parameters: 

Cartesian  Grid  Module  -EH 

a.  Save  Bathymetry  Contours  in  a  DXF  file: 


1.  Edit  Window  |  Solution  ->  Generic  Data  sets. 

2.  Display  |  Display  Options 
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3.  Check,  Contours  ->  Options  ->  Contour  Method  ->  Normal 
linear  Contours  ->  OK. 

4.  OK. 

5.  File  |  SaveAs  ->  Bathymetry  File  ->  File  Type  =  .DXF 

b.  Load  DXF  bathymetry  contours  in  SMS: 

1.  File  |  Open  ->  Bathymetry  File 

2.  Open  File  Format  dialog  ->  Open  File  As  ->  .DXF,  OK. 

3.  Dialog,  DXF  Data  Exists.  Replace  or  Append?  Replace. 
OK. 

c.  Display  bathymetry  over  model  wave  estimates: 

1.  Edit  Window  |  Solution  ->  project(STWAVE) 

2.  Display  |  DisplayOptions 

3.  Check,  Contours  ->  Options  ->  Contour  Method  ->Color 
Filled  between  Contours  ->  OK. 

4.  OK. 

Display  monitoring  stations 

Menu  Bar 

a.  Display  |  Display  Options 

b.  Check,  Cells 

c.  Click,  Cell  Options  ->  Check,  Observation  Symbol  ->  OK. 

d.  OK. 

Visualize  calculated  wave  height,  period,  or  direction  fields 

Cartesian  Grid  Module  -H 

a.  Select  an  active  data  set  by  setting  these  three  lists  in  the  edit  window: 
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1. 


Scalar 


2.  Vector 

3.  Time-Step 


Visualize  calculated  wave  height,  period,  or  direction  cell 

Follow  previous  instructions  for  “Visualize  calculated  wave  height,  period, 
or  direction  fields,”  then: 

Select,  Cell  Tool  -  B. 

a.  Click  the  mouse  cursor  anywhere  on  the  display  of  the  calculated  wave 

field. 

1 .  The  cell  identifier  (I,J)  coordinate  and  other  information  for  the  cell 
is  given  at  the  bottom  left  portion  of  the  SMS  screen. 

2.  The  scalar  value  for  the  cell  of  the  active  data  set  is  displayed  in  the 
edit  window  after  “S:”.  “S:”  is  located  directly  below  “Scalar”  in  the 
edit  window. 

3.  Because  this  simulation  includes  currents,  the  active  data  set  values 
for  “Vx”  and  “Vy”  are  displayed  to  the  right  of  “S:”. 


Dynamically  visualize  calculated  wave  height,  period,  or  direction 
cell 


Follow  previous  instructions  for  “Visualize  calculated  wave  height,  period, 
or  direction  fields,”  then: 

Select,  Cell  Tool  IH. 


a.  As  the  mouse  cursor  is  moved  anywhere  over  the  estimated  wave  field 
display,  notice  the  third  line  of  the  edit  window  displaying  x-,  y-,  z- 
coordinates  and  model  results  corresponding  to  the  active  data  set. 


Visualize  observation  spectra  at  monitoring  station 

Energy  spectra  saved  at  monitoring  stations  can  be  visualized  from  the 
Spectral  Energy  Dialog  (Figure  9).  Use  the  Spectral  Energy  dialog’s  edit  window 

[Xy  1 

or  data  browser  IHl  to  select  a  monitoring  station  and  time-step.  Monitoring 
stations  are  identified  by  their  I,J  grid  location  (e.g.,  Node  l  10157, 

Node_131_l  17,  Node_146_86).  Polar  and  Cartesian  contour  options  and  energy 
frequency  and  direction  plots  are  available  for  displaying  spectra. 
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Spectral  Energy  Dialog 


View: 

O  Polar 
Cartesian 

Selected  Spectrum: 
Wind  Speed  (m/s): 


Wind  Dir.  (deg): 


10.000 

Tidal  Elev.  (m): 


r  Freq.  Plot 
r  Dir.  Plot 


Figure  9.  Spectral  energy  dialog 

Cartesian  Grid  Module  -  S 

a.  Stwave  |  Spectral  Energy 

b.  Function:  Node_110_157 

c.  Time-Step:  200700 

d.  Under  the  heading.  View,  check:  Cartesian 

e.  Under  the  heading,  Display ,  check:  Freq.  Plot 
f  Under  the  heading,  Display,  check:  Dir.  Plot 
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Animation  film  loops 

Multiple  STWAVE  simulations  produce  multiple  solutions.  These  solutions  can 
be  animated  and  saved  as  an  AVI  file.  Refer  to  the  SMS  “Animation”  Help  topic 
for  additional  information. 

Cartesian  Grid  Module  -H 

a.  Data  |  Film  Loop... 

b.  Setup 

c.  On  the  Film  Loop  Options  dialog  in  the  Data  Options  section,  select: 

1.  Scalar  data  set 

2.  Run  simulation  from  time:  175500 

3.  Run  simulation  to  time:  200700 

4.  Match  time-steps 

5.  OK 

Two-Dimensional  Plots 

SMS  offers  various  2-D  plots.  Refer  to  the  SMS  Help  topic  on  2-D  plots  for  a 
complete  list. 

A  profile  plot  displays  the  variation  of  one  or  more  scalar  data  sets  associated 
with  a  mesh  or  grid  along  arcs  in  the  observation  coverage.  For  instance,  wave 
heights  along  an  arc  may  be  displayed  for  a  specific  event  identifier  (Figure  10). 
The  arc  can  be  drawn  anywhere  on  the  data  set: 

Map  Module  -ES 

a.  Feature  Objects  |  Coverages... 

1.  Coverage  Type ->  Observation 

b.  Display  |  Plot  Options 

1.  New 

2.  Plot  Type:  Observation  Profile 
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3.  Options: 

(a)  Select,  Use  selected  data  sets: 


(1)  Highlight:  Height 

(2)  Use  time-step:  200700 

(3)  OK 

4.  Labels,  scales... 

(a)  Title  2:  Event  Identifier:  200700  -  Units:  meters 

(b)  X-axis:  Distance 

(c)  Y-axis:  Wave  Height 

(d)  OK 


5.  OK 


c.  Create  Feature  Arc  Tool 


I .  Click  the  mouse  on  the  data  set  where  you  want  to  begin  the  arc. 
Move  and  click  the  mouse  to  draw  an  arc.  Double-click  the  mouse  to 
end  the  arc. 


d.  Display  |  Plot  Window 

1 .  A  plot  window  opens  below  the  edit  window  and  above  the 

ST  WAVE  data  set  window.  You  can  enlarge  the  plot  window  by 
placing  the  mouse  cursor  at  the  bottom  of  the  window.  When  the 
cursor  changes  to  a  bar  with  an  arrow-  head  on  each  end,  press  the 
mouse  button  and  drag  the  window  to  the  desired  size. 
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6  Example  Applications 


In  this  chapter,  three  example  applications  of  STWAVE  are  given.  The 
purpose  of  these  examples  is  to  demonstrate  how  the  model  can  be  applied  and 
provide  sample  problems  that  a  new  user  can  follow  to  learn  the  mechanics  of 
model  application.  The  first  example  is  wave  transformation  and  wave-current 
interaction  at  a  complex  coastal  entrance  on  the  East  Coast.  The  second  example 
is  wind-wave  generation  over  a  simple  bathymetry.  The  third  example  is  wave 
transformation  at  a  West  Coast  entrance  channel  using  the  files  generated  in  SMS 
following  the  procedure  outlined  in  Chapter  5.  The  first  and  third  examples 
focus  on  wave  propagation  and  the  second  example  on  wave  generation  (using 
the  source  term  options). 


Example  1 :  Wave  Propagation  at  Ponce  de  Leon 
Inlet 

The  first  example  is  based  on  a  field  site,  Ponce  de  Leon  Inlet,  located  on  the 
east  coast  of  Florida.  The  nearshore  bathymetry  at  Ponce  Inlet  is  complex 
because  of  the  presence  of  the  inlet,  a  navigation  channel,  jetties,  and  a  large  ebb 
shoal.  Tidal  currents  at  the  entrance  can  exceed  1  m/sec.  Engineering  problems 
at  Ponce  Inlet  that  require  information  about  the  nearshore  wave  field  include 
scour  in  the  inlet  throat  near  the  north  jetty,  erosion  of  the  north  spit  in  the 
interior  of  the  inlet,  migration  of  the  navigation  channel,  and  possible  breaching 
shoreward  of  the  north  jetty.  More  information  about  Ponce  Inlet  is  given  in  a 
reconnaissance  report  by  the  U.S.  Army  Engineer  District,  Jacksonville  (1993). 
Also,  a  field-data  collection  program  conducted  at  Ponce  Inlet  from  September 
1995  through  September  1997  is  documented  by  King  et  al.  (1999).  STWAVE 
applications  to  Ponce  Inlet  are  documented  in  Smith,  Militello,  and  Smith  (1998) 
and  Smith  and  Smith  (in  preparation). 


Model  input 

Model  parameters.  Figure  1 1  shows  the  bathymetry  contours  within  the 
STWAVE  modeling  domain  selected  for  Ponce  Inlet.  The  offshore  boundary 
location  is  selected  at  an  approximate  depth  of  1 8  m.  This  depth  is  chosen 
because  the  bottom  contours  are  fairly  straight  and  parallel  at  this  depth  (out  of 
the  influence  of  the  ebb  shoal)  and  because  a  wave  gauge  (designated  as  DWG1 
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in  Figure  11)  was  located  near 
this  contour  to  provide  the  input 
wave  conditions.  The  lateral 
boundaries  of  the  domain  are 
positioned  away  from  the 
influence  of  the  inlet,  to  areas  of 
fairly  straight  and  parallel  depth 
contours.  A  50-m  grid  cell 
spacing  is  selected  to  resolve  the 
ebb  shoal  and  inlet  bathymetry. 
To  cover  the  domain  with  50-m 
resolution  requires  53  cells 
across  the  shore  (NI)  and 
112  cells  alongshore  (NJ). 
Typical  wave  periods  for  this 
site  are  15  to  5  sec  (0.0667  to 
0.2  Hz).  To  resolve  this  range, 
30  frequency  bins  are  used  with 
an  initial  frequency  of  0.03 1  Hz 
and  a  frequency  increment  of 
0.0078  Hz  (range  of  frequencies 
is  0.031  to  0.258  Hz).  Because 
the  wave  propagation  distances 
on  the  Ponce  Inlet  grid  are  short 
(2  km),  source  terms  are 
neglected  (IPRP  =1).  Even 
with  strong  winds,  wave  growth 
over  the  short  distance  would  be 
small.  Wave-current  interaction, 
though,  could  be  significant  near 
the  inlet,  so  ICUR  =  1  to  include 
wave-current  interaction  is 
selected.  Breaker  indices  are  not 
saved  (IBREAK  =  0),  and 
neither  are  radiation  stress 
gradients  (IRS  =  0).  Selected 
model  output  locations  are 
chosen  to  match  three  directional 


Figure  1 1 .  Ponce  Inlet  bathymetry  and 
gauge  positions 


wave  gauge  positions  at  Ponce 

Inlet  (see  Figure  1 1).  The  depths  at  the  three  gauges  are  approximately  14,  7,  and 
5  m.  The  gauges  are  located  near  the  offshore  boundary,  on  the  ebb  shoal,  and  in 
the  outer  inlet  throat  (north  of  the  ebb  shoal). 


Bathymetry.  The  bathymetry  contours  for  Ponce  Inlet  are  shown  in 
Figure  1 1 .  The  grid  was  developed  using  the  SMS  grid  generator,  which  accepts 
random  x,  y,  and  depth  triplets  and  interpolates  them  onto  a  Cartesian  grid  with  a 
given  origin,  orientation,  and  resolution.  In  this  case,  high-resolution  bathymetry 
near  the  inlet  was  obtained  from  a  Scanning  Hydrographic  Operational  Airborne 
Lidar  System  (SHOALS)  (Lillycrop,  Parson,  and  Irish  1996)  survey  and 
supplemented  with  National  Oceanic  and  Atmospheric  Administration  (NOAA) 
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digital  bathymetry  (used  to  define  offshore  and  bay  bathymetry)  and  shorelines. 
The  north  jetty  at  Ponce  Inlet  is  represented  in  the  model  as  a  series  of  land  cells. 
The  jetty  is  specified  as  being  three  grid  cells  wide,  which  is  wider  than  the  real 
structure.  Three  cells  are  required  to  block  wave  energy  from  propagating 
through  the  structure.  Future  upgrades  to  the  model  will  provide  the  capability  to 
represent  structure  widths  that  are  less  than  the  grid  cell  spacing. 

Incident  wave  spectrum.  The  incident  wave  spectrum  for  Ponce  Inlet  was 
generated  using  a  TMA  spectral  shape  (with  a  spectral  peakedness  parameter,  y- 
5.0),  cos™ a  directional  distribution  (with  nn  =  12),  Hm0  =  5.2  m,Tp-  12.8  sec, 
and  am  =  -7  deg  (see  Table  1  for  guidance  on  spectra  shape  and  directional 
spreading  parameters).  These  incident  conditions  represent  extreme  high  wave 
conditions  measured  at  Ponce  Inlet  in  March  1 996. 

Current  field.  The  current  field  for  Ponce  Inlet  was  generated  using  the 
tidal  circulation  model  ADCIRC  (Luettich,  Westerink,  and  Scheffner  1992).  The 
current  field  output  from  ADCIRC  was  interpolated  onto  the  STWAVE  grid 
using  SMS.  The  current  field  is  plotted  in  Figure  12. 


Results 

The  wave  parameters  for  the  three  gauge  positions  at  Ponce  Inlet  are  given  in 
the  selhts  file: 


Date 

I 

J 

77mo(m) 

Tp{  sec) 

am (deg) 

96031202 

15 

21 

5.25 

12.8 

-6. 

96031202 

15 

65 

4.13 

12.8 

-5. 

96031202 

28 

49 

2.86 

12.8 

18. 

At  the  first  gauge  position  (I  =  15,  J  =  21),  the  wave  height  has  shoaled  slightly 
(from  5.2  m  at  the  boundary  to  5.25  m)  and  turned  slightly  more  shore  normal 
(from  -7  deg  at  the  boundary  to  -6  deg).  At  the  second  gauge  (I  =  15,  J  =  65), 
located  on  the  ebb  shoal,  the  wave  energy  has  dissipated  because  of  breaking 
(reduction  of  21  percent  in  wave  height  from  the  offshore  boundary).  At  the 
most  shoreward  gauge  (I  =  28,  J  =  49),  the  energy  has  been  dissipated 
significantly  (45  percent  reduction  in  wave  height  from  the  offshore  boundary) 
and  the  mean  direction  has  refracted  from  -7  deg  at  the  boundary  to  1 8  deg  in  the 
outer  throat.  The  wave  directions  at  this  gauge  are  turning  to  align  normal  to  the 
ebb  shoal  that  is  located  to  the  south  of  the  gauge  (positive  angles  are  more 
southerly  directed  and  negative  angles  are  more  northerly  directed). 

The  trends  of  wave-height  reduction  and  turning  of  the  wave  angle  can  also 
be  illustrated  by  examining  the  one-dimensional  wave  spectra  plotted  in 
Figure  13.  The  shape  of  the  spectra  stays  quite  similar,  but  the  energy  is  reduced 
because  of  depth-limited  breaking  from  the  offshore  to  the  outer  throat  gauge. 
Another  way  to  view  the  spectrum  is  to  integrate  over  all  frequencies  to  examine 
the  directional  distribution  of  the  wave  energy  (Figure  14).  The  directional 
distribution  narrows  as  well  as  reduces  in  energy  between  the  offshore  and  ebb 
shoal  gauges  because  of  refraction.  Between  the  offshore/ebb  shoal  gauges  and 
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the  outer  throat  gauge,  the  mean  direction  shifts  by  25  deg,  again  because  of 
refraction. 

Figure  15  shows  a  contour  plot  of  the  wave  heights  over  the  entire  STWAVE 
domain.  The  gray  shade  contours  represent  wave  height,  and  the  background  line 
contours  represent  the  bathymetry.  For  these  incident  wave  conditions,  the 
wave-height  variations  closely  follow  the  depth  contours,  because  depth-limited 
breaking  is  the  dominant  process  and  thus  the  wave  height  contours  mimic  the 
depth  contours.  For  cases  with  less  extreme  incident  wave  heights,  the  maximum 
wave  heights  are  found  on  the  ebb  shoal,  where  energy  is  focused  by  refraction. 
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Example  2:  Wind-Wave  Generation 

For  the  second  example,  consider  a  bay  that  is  rectangular  in  shape,  20  km 
east-west  and  10  km  north-south.  For  practical  purposes,  the  bay  is  isolated  from 
ocean  waves.  The  water  depth  is  a  constant  15  m,  with  a  gentle  side  slope  of 
1 : 1 00  at  the  eastern  end  of  the  bay  and  vertical  walls  on  the  other  three  sides. 

The  bay  has  a  rectangle  shape  with  simplified  bathymetry,  so  users  can  easily 
develop  the  input  bathymetry.  The  grid  domain  is  shown  in  Figure  16.  The  area 
of  interest  is  the  eastern  shoreline  of  the  bay,  where  a  project  is  planned  to  mine 
sand  to  construct  a  beach.  For  a  project  such  as  this,  wave  information  may  be 
required  to  evaluate  sediment-transport  rates  at  the  borrow  site  to  estimate 
infilling  rates  or  evaluate  longshore  and  cross-shore  sediment-transport  rates  at 
the  newly  constructed  beach  to  estimate  the  project’s  life  and  renourishment 
requirements. 


Model  input 

Model  parameters.  Both  wave  propagation  and  source  terms  (IPRP  =  0)  are 
required  for  this  example,  but  currents  are  negligible  (ICUR  =  0).  Breaker 
indices  and  radiation  stress  gradients  are  not  required  (IBREAK  =  0  and  IRS  = 
0).  The  special  output  points  are  selected  to  examine  wave  growth  along  the 
major  axis  of  the  bay.  The  borrow  site  is  located  near  grid  cell  (85,25).  The 
model  parameter  input  file  is: 

0  0  0  0  10 
5  25 
25  25 
50  25 
75  25 
85  25 
95  25 

97  25 

98  25 

99  25 

100  25 


Chapter  6  Example  Applications 


55 


Bathymetry.  The  bathymetry  and  shoreline  are  very  simple,  and  fine 
bathymetric  resolution  is  not  required.  A  grid  spacing  (DXINC)  of  200  m  is 
selected,  with  101  cells  in  the  cross-shore  (NI)  and  51  cells  in  the  alongshore 
(NJ).  To  simplify  the  input  file  for  this  example,  the  bay  is  assumed  to  have 
vertical  walls  on  its  north,  south,  and  west  shores  and  a  1:100  slope  on  the 
eastern  shore.  The  depth  in  all  grid  cells  along  the  entire  grid  boundary  is  set  to  - 
1  for  land.  The  first  two  grid  rows  (J  =  5 1  and  J  =  50)  of  the  depth  file  are  given 
as  follows: 
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The  full  depth  file  would  repeat  the  second  series  for  rows  J  =  49  through  2,  and 
the  depths  for  row  J  =  1  would  again  all  be  -1 .,  identical  to  the  first  series  in  the 
file.  Note  that  the  first  line  in  the  depth  file  represents  the  northern  edge  of  the 
grid  (J  =  NJ),  reading  from  west  to  east  (I  =  1,NT),  and  the  last  series  in  the  files 
would  be  the  southern  edge  of  the  grid  (J  =  1). 

Incident  wave  spectrum.  Twenty  frequency  bins  are  used  with  an  initial 
frequency  of  0.04  Hz  and  a  frequency  increment  of  0.04  Hz.  This  coarse 
frequency  resolution  is  sufficient  because,  even  at  high  wind  speeds,  the  20-km 
fetch  will  limit  peak  wave  periods  to  approximately  5  to  6  sec.  The  input 
spectrum  is  very  simple  for  this  example  case.  Because  the  bay  is  treated  as  an 
enclosed  basin,  the  spectrum  on  the  “offshore”  boundary  (the  western  land 
boundary  for  this  case)  is  zero.  The  input  spectrum  includes  the  header  with  the 
date  and  wind  information,  followed  by  700  zeros  (20  frequencies  by  35 
directions).  For  a  date  identifier  of  980923,  a  wind  speed  of  20  m/sec,  a  wind 
direction  of  0  deg  relative  to  the  STWAVE  grid  (wind  blowing  from  the  west,  for 
this  example),  a  peak  frequency  of  0.8  Hz  (set  to  the  highest  frequency  because 
there  is  no  energy  in  the  input  spectrum),  and  0.0  m  water  elevation  correction; 
the  beginning  of  the  spectrum  file  is  given  by: 

20  35 

0.04  0.08  0.12  0.16  0.20  0.24  0.28  0.32  0.36  0.40 
0.44  0.48  0.52  0.56  0.60  0.64  0.68  0.72  0.76  0.80 
980923  20.  0.  0.8  0. 

In  the  input  file,  these  lines  would  be  followed  by  700  zeros  (35  directions  by 
20  frequencies)  to  represent  an  incident  spectrum  with  no  energy. 
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Current  field.  No  current  field  file  is  specified  for  this  case. 


Results 

The  results  of  greatest  interest  for  this  case  are  the  wave  parameters  at  the 
proposed  borrow  site  (I  =  85,  J  =  25).  The  wave  height  is  1 .52  m,  the  peak 
period  is  4.6  sec,  and  the  direction  is  0  deg  (from  the  west,  for  this  grid 
orientation).  The  output  file  with  selected  wave  parameters  is  given  as  follows 
(header  has  been  added): 


Date 

I 

J 

Tp{  sec) 

ctm(deg) 

980923 

25 

25 

0.86 

3.2 

0. 

980923 

50 

25 

1.22 

3.9 

0. 
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75 

25 

1.45 

4.4 

0. 
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85 

25 

1.52 

4.6 

0. 

980923 

95 

25 

1.55 

4.7 

0. 

980923 

96 

25 

1.53 

4.7 

0. 

980923 

97 

25 

1.50 

4.7 

0. 

980923 

98 

25 

1.47 

4.8 

0. 

980923 

99 

25 

1.23 

4.8 

0. 

980923 

100 

25 

0.59 

4.8 

0. 

These  output  wave  parameters  across  the  long  axis  of  the  bay  demonstrate 
the  wave  growth  as  a  function  of  distance  along  the  fetch.  These  wave  heights 
and  periods  are  plotted  as  a  function  of  distance  in  Figure  17.  The  wave  height 
grows  with  fetch  until  the  depth  decreases  on  the  eastern  shore  of  the  bay,  and  the 
wave  height  decays  because  of  depth-limited  breaking.  The  wave  direction  down 
the  long  axis  of  the  bay  is  aligned  with  the  wind  direction  (0  deg).  Figure  18 
shows  selected  one-dimensional  spectra  along  the  long  axes  of  the  bay  (legend 
provides  I,  J  locations).  Moving  from  west  to  east  (I  =  25  to  95),  the  peak 
frequencies  decrease  (peak  periods  increase)  and  the  total  energy  increases,  as  the 
waves  grow  from  input  wind  energy.  Between  I  =  95  and  100,  the  energy  density 
decreases  because  of  depth-limited  breaking.  These  one-dimensional  spectra 
were  calculated  from  the  spec .  out  file  by  integrating  the  spectra  over  all 
directions. 

To  evaluate  the  impacts  of  the  borrow  site  on  the  shoreline  or  the  sediment 
infilling  rate  at  the  borrow  site,  additional  STWAVE  runs  would  be  required. 
Typically,  the  local  wave  climate  would  be  examined,  and  representative  as  well 
as  maximum  wind  conditions  would  be  used  as  input  for  multiple  model  runs. 
Results  from  these  runs  would  then  be  used  statistically  to  evaluate  sediment- 
transport  potential  at  the  borrow  site  and  along  the  shoreline.  Note  that  the 
strongest  wind  speed  does  not  necessarily  lead  to  the  largest  waves  or  largest 
sediment-transport  rates.  Wind  direction  is  also  a  critical  parameter. 

This  example  of  wave  growth  was  for  a  completely  enclosed  basin.  The 
wave  growth  capabilities  of  the  model  are  also  applicable  on  the  open  coast.  For 
open-coast  applications,  the  input  spectrum  would  likely  not  equal  zero,  and  the 
offshore  and  lateral  grid  boundaries  would  not  be  land.  For  open  coastal 


Chapter  6  Example  Applications 


57 


Figure  17.  Wave  height  and  period  as  a  function  of  distance  for  20  m/sec  wind 
speed 


Figure  18.  One-dimensional  wave  spectra  for  20  m/sec  wind  speed 
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Peak  Wave  Period  (sec) 


applications,  the  local  wave  generation  (represented  by  the  source  terms)  is  often 
an  essential  process,  increasing  the  wave  height  near  shore  and  altering  wave 
directions. 


Example  3:  Wave  Propagation  at  Yaquina  Bay 
Entrance,  OR 

The  final  example  is  based  on  a  field  site,  Yaquina  Bay  Entrance,  OR.  The 
nearshore  bathymetry  at  Yaquina  is  complex  because  of  the  presence  of  the 
entrance  channel,  rock  reefs,  and  jetties.  Tidal  currents  at  the  entrance  exceed 
1 .5  m/sec,  and  the  tidal  range  is  greater  than  2  m.  The  wave  climate  off  the 
Oregon  coast  is  severe  with  observed  deepwater  wave  heights  exceeding  12  m. 
The  north  Yaquina  jetty  was  repaired  in  1988  and  has  experienced  progressive 
loss  of  the  head  of  the  jetty  since  1996  (approximately  110  m  of  jetty  length  has 
been  lost).  Quantitative  information  about  waves  at  the  entrance  was  required  to 
support  structure  redesign  and  to  evaluate  navigability.  More  information  about 
the  Yaquina  Bay  Entrance  north  jetty  is  given  by  Hughes  et  al.  (1995). 

The  nearest  wave  gauge  to  Y aquina  Bay  is  a  National  Data  Buoy  Center 
buoy  located  35  km  off  the  coast  in  a  depth  of  130  m.  It  is  not  practical  to 
transform  the  measured  waves  from  the  buoy  to  the  entrance  at  fine  resolution. 
Although  high  resolution  (25  m)  is  required  in  the  nearshore  to  define  the  reefs 
and  entrance  channel,  coarser  resolution  is  sufficient  offshore  (where  high- 
resolution  bathymetry  is  not  available).  For  this  application,  a  fine  nearshore  grid 
was  nested  with  a  coarse  offshore  grid.  The  offshore  grid  has  a  resolution  of 
250  m  with  125  cells  across  the  shore  and  253  cells  along  the  shore.  The  grid 
extends  from  an  offshore  depth  of  130  m  (at  the  buoy  location)  to  a  nearshore 
depth  of  47  m.  NOAA  bathymetry  was  used.  Both  the  offshore  and  nearshore 
grids  were  aligned  with  the  nearshore  contours  (x-axis  at  350  deg  counterclock¬ 
wise  from  west),  so  no  rotation  of  the  spectra  was  required  between  offshore  and 
nearshore  simulations.  Measured  heights,  periods,  and  directions  from  the 
offshore  buoy  were  then  transformed  across  the  offshore  grid  (including  wind 
input,  but  excluding  currents  and  tide  variations).  Output  spectra  were  saved  in 
the  most  landward  column  and  used  as  the  input  spectra  for  the  nearshore  grid. 
The  nearshore  model  runs  are  described  in  the  following  paragraphs. 


Model  input 

Model  parameters.  Figure  19  shows  the  bathymetry  contours  within  the 
STWAVE  nearshore  modeling  domain  selected  for  the  Yaquina  Bay  Entrance. 
The  offshore  boundary  location  is  selected  at  an  approximate  depth  of  47  m. 

This  depth  is  chosen  because  the  bottom  contours  are  fairly  straight  and  parallel 
at  this  depth  (out  of  the  influence  of  the  ebb  shoal).  The  lateral  boundaries  of  the 
domain  are  positioned  away  from  the  influence  of  the  entrance,  to  areas  of  fairly 
straight  and  parallel  depth  contours.  A  25-m  grid  cell  spacing  is  selected  to 
resolve  the  reefs  and  entrance  bathymetry.  To  cover  the  domain  with  25-m 
resolution  requires  188  cells  across  the  shore  (NI)  and  236  cells  alongshore  (NJ). 
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Typical  wave  periods  for  this  site  are  20  to  6  sec  (0.05  to  0.167  Hz).  To  resolve 
this  range,  30  frequency  bins  are  used  with  an  initial  frequency  of  0.02  Hz  and  a 
frequency  increment  of  0.01  Hz  (range  of  frequencies  is  0.02  to  0.3 1  Hz).  The 
same  frequency  range  and  resolution  were  used  in  both  the  offshore  and  near¬ 
shore  simulations,  so  no  interpolation  was  required  between  offshore  grid  output 
spectra  and  the  nearshore  grid  input  spectra.  Because  the  wave  propagation 
distances  on  the  nearshore  grid  are  relatively  short  (5  km),  source  terms  are 
neglected  (fPRP  =1).  Wave-current  interaction,  though,  is  significant  near  the 
inlet,  so  ICUR  =  1  to  include  wave-current  interaction  was  selected.  Breaker 
indices  and  radiation  stress  gradients  are  saved  (IBREAK  =  1  and  IRS  =  1). 

Bathymetry.  The  bathymetry  contours  for  Yaquina  Bay  are  shown  in 
Figure  19.  The  grid  was  developed  using  the  SMS  grid  generator  as  presented  in 
Chapter  5.  The  Portland  District  provided  high-resolution  bathymetry  near  the 
entrance  from  multibeam  fathometry.  These  data  were  supplemented  with 
NOAA  bathymetry  and  shorelines.  The  Yaquina  jetties  are  represented  in  the 
model  as  a  series  of  land  cells.  The  jetty  is  specified  as  being  three  grid  cells 
wide,  which  is  wider  than  the  real  structure.  Three  cells  are  required  to  block 
wave  energy  from  propagating  through  the  structure. 


Figure  19.  Nearshore  Yaquina  Bay  Entrance  bathymetry  (from  SMS) 
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Incident  wave  spectrum.  The  incident  wave  spectrum  for  Yaquina  Bay  was 
generated  using  a  TMA  spectral  shape  (with  a  spectral  peakedness  parameter, 
y=  8),  cos’" a  directional  distribution  (with  nn  =  30),  Hm0  =  12.5  m,  Tp  =  20  sec, 
and  am  =  55  deg  (from  the  Southwest)  on  the  offshore  grid.  The  nearshore  grid 
was  forced  with  the  output  spectrum  at  a  depth  of  47  m  from  the  offshore  grid. 
These  incident  conditions  represent  extreme  high  wave  conditions  from  a  storm 
in  1999. 

Current  field.  The  current  field  for  Yaquina  Bay  was  generated  using  the 
tidal  circulation  model  ADCIRC.  The  current  field  output  from  ADCIRC  was 
interpolated  onto  the  STWAVE  grid  using  SMS  as  described  in  Chapter  5.  The 
ebb  and  flood  current  fields  are  plotted  in  Figures  17  and  18  for  the  area  near  the 
jetty  tips.  The  current  fields  show  not  only  the  tidal  current  through  the  entrance, 
but  also  a  longshore  current  toward  the  north  generated  by  strong  alongshore 
storm  winds.  Representative  nearshore  tide  levels  of -0.90  m  (at  peak  ebb)  and 
0.61  m  (at  peak  flood)  from  ADCIRC  were  applied  in  STWAVE. 


Results 

The  wave  height  contours  in  Figures  20  and  21  show  that  the  waves  in  the 
nearshore  in  this  extreme  storm  event  are  breaking  due  to  depth  limitations  (note 
gradient  in  the  wave  height  contours).  The  1.5-m  higher  water  level  on  flood 
accounts  for  the  larger  wave  heights  during  flood  than  during  ebb  (generally, 
larger  wave  heights  are  expected  on  ebb  due  to  wave-current  interaction,  unless 
the  wave  heights  are  depth  limited,  as  they  are  in  this  example).  Figure  22 
(generated  in  SMS)  shows  the  difference  in  wave  height  (flood  minus  ebb)  for 
this  example.  The  red/orange  contours  show  areas  where  the  wave  height  is  1 .5 
to  2  m  higher  on  flood  (at  the  jetty  tips).  The  blue  contours  show  areas  where  the 
wave  height  is  0.5  to  1 .0  m  higher  on  ebb  (in  the  interior  channel,  where  the 
water  depth  is  greater  and  the  current  stronger). 
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Figure  20.  Ebb  current  vectors  and  wave  height  contours  near  the  Yaquina  jetty 
tips  (incident  conditions:  Hmo  =  12.5  m,  Tp  =  20  sec,  waves  from  the 
southwest,  tide  =  -0.90  m  mean  tide  level) 


Figure  21 .  Flood  current  vectors  and  wave  height  contours  near  the  Yaquina 

jetty  tips  (incident  conditions:  Hmo  =  12.5  m,  Tp  =  20  sec,  waves  from 
the  southwest,  tide  =  0.61  m  mean  tide  level) 
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Figure  22.  Wave  height  difference  (flood  -  ebb)  near  the  Yaquina  jetty  tips 
(incident  conditions:  Hmo  =  12.5  m,  Tp  =  20  sec,  waves  from  the 
southwest) 
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Appendix  A 

Sample  Model  Parameter  File 


The  parameters  required  to  run  STWAVE  are  specified  in  the 
options  .  std  file.  A  description  of  the  file  and  formatting  is  given  in  Chapter 
3  of  the  main  text.  An  example  parameter  file  for  Yaquina  Bay  Entrance  at 
Newport,  OR  (see  Example  1  in  Chapter  6  -  main  text)  is  given  as  follows: 

0  1113 

110  157 
131  117 
146  86 


Appendix  A:  Sample  Model  Parameter  File 


A1 


Appendix  B 

Sample  Bathymetry  File 


The  bathymetry  required  to  run  STWAVE  is  specified  in  the  depth .  in 
file.  A  description  of  the  file  and  formatting  is  given  in  Chapter  3, 

“Bathymetry.”  An  example  of  the  first  row  (starting  with  J  =  NJ)  of  the 
bathymetry  file  for  Yaquina  Bay  Entrance  at  Newport,  OR  (Chapter  5,  “Build 
STWAVE  Cartesian  Grid”)  is  given  as  follows  (depths  in  meters  relative  to  mean 
tide  level): 


188  236  25.000000 


46.310997 

46.545502 

46.476299 

46.465614 

46.099133 

45.525253 

45.168720 

44.934937 

44.753963 

45.227501 

45.381058 

45.078934 

44.749165 

44.173389 

43.706909 

43.200413 

42.864304 

42.718304 

42.552567 

41.984577 

41.584068 

41.581631 

41.595505 

41.608906 

41.114693 

40.779945 

40.480919 

39.474430 

39.658737 

39.700066 

39.478432 

39.481739 

39.362648 

38.426922 

37.489922 

36.998165 

36.513798 

36.748604 

37.013935 

36.346466 

35.320770 

35.138374 

35.586994 

35.580254 

35.479980 

34.844318 

34.245529 

33.885906 

33.983131 

34.072853 

33.580402 

33.224022 

32.847286 

32.452747 

32.628101 

32.529854 

31.978901 

31.437374 

30.874306 

30.525604 

30.289108 

30.109001 

29.668436 

29.209021 

29.025633 

28.913048 

28.339758 

27.856483 

27.660454 

27.450241 

27.198229 

26.932903 

26.698179 

26.465672 

26.488281 

26.445293 

25.880035 

25.253780 

24.910845 

24.707207 

24.630911 

24.686291 

24.185822 

23.762606 

23.608511 

23.580662 

23.160681 

22.718637 

22.439632 

22.313131 

21.931955 

21.434870 

21.035362 

20.962799 

20.753426 

20.374048 

20.099985 

20.066952 

20.025595 

19.114258 

18.769199 

19.031534 

18.911980 

17.971195 

17.398960 

17.377380 

17.077021 

16.958942 

16.618668 

16.451231 

16.231750 

15.705775 

15.137889 

14.560424 

14.812372 

15.064322 

15.316270 

13.664304 

16.208809 

16.184008 

16.021345 

15.079105 

14.334208 

14.016557 

13.685314 

13.354072 

13.022829 

12.641250 

12.542627 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

-2.000000 

-2.000000 

-2.000000 

-2.000000 

-2.000000 

-2.000000 

-2.000000 

-2.000000 

-2.000000 

-2.000000 

-2.000000 

-2.000000 

-2.000000 

-2.000000 

-2.000000 

Appendix  B:  Sample  Bathymetry  File 
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Appendix  C 
Sample  Incident  Wave 
Spectrum  File 


The  input  two-dimensional  spectrum  required  to  run  STWAVE  is  specified 
in  the  spec .  in  file.  A  description  of  the  file  and  formatting  is  given  in 
Chapter  3,  “Incident  Wave  Spectra.”  An  example  of  a  file  header  and  the  energy 
densities  (meters  squared/hertz/radian)  in  the  first  five  frequencies  (0.01000  to 
0.05000  Hz,  with  35  direction  bins  per  frequency)  of  the  spectrum  for  Yaquina 
Bay  Entrance  at  Newport,  OR  (Chapter  5,  “Generate  Input  Spectra  for 
STWAVE”)  are  given  as  follows: 

25  35 

0.01000  0.02000  0.03000  0.04000  0.05000  0.06000  0.07000  0.08000  0.09000  0.10000 
0.11000  0.12000  0.13000  0.14000  0.15000  0.16000  0.17000  0.18000  0.19000  0.20000 


0.21000  0.22000  0.23000  0.24000  0.25000 


200700  15 

.000000  10 

.000000  0, 

.070000  -1, 

.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000005 

0.000035 

0.000195 

0.000824 

0.002753 

0.007495 

0.016990 

0.032582 

0.053484 

0.075781 

0.093203 

0.099824 

0.093203 

0.075781 

0.053484 

0.032582 

0.016990 

0.007495 

0.002753 

0.000824 

0.000195 

0.000035 

0.000005 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000003 

0.000058 

0.000688 

0.005347 

0.029759 

0.125738 

0.420199 

1.144050 

2.593323 

4.973229 

8.163601 

11.566981 

14.226231 

15.236800 

14.226231 

11.566982 

8.163602 

4.973230 

2.593323 

1.144051 

0.420199 

0.125738 

0.029759 

0.005347 

0.000688 

Appendix  C:  Sample  Incident  Wave  Spectrum  File 


Appendix  D 

Sample  Current  Field  File 


The  input  current  field  required  to  run  STWAVE  is  specified  in  the 
current .  in  file.  A  description  of  the  file  and  formatting  is  given  in 
Chapter  3,  “Incident  Wave  Spectra.”  An  example  of  the  first  row  (J  =  NJ)  of  the 
current  file  for  Yaquina  Bay  Entrance  at  Newport,  OR  (Chapter  5,  “Wave- 
Current  Interaction”)  is  given  as  follows  ( u  and  v  velocity  components  in 
meters/second): 


188  236  25.000000 
200700 


-0.074441 

0.521558 

-0.075095 

0.522148 

-0.075748 

0.522737 

-0.076402 

0.523327 

-0.077055 

0.523917 

-0.077708 

0.524506 

-0.078377 

0.525387 

-0.079053 

0.526411 

-0.079728 

0.527434 

-0.080404 

0.528457 

-0.081080 

0.529480 

-0.081755 

0.530503 

-0.080008 

0.529364 

-0.079146 

0.528634 

-0.078284 

0.527905 

-0.077422 

0.527175 

-0.076560 

0.526446 

-0.075699 

0.525716 

-0.074837 

0.524986 

-0.073975 

0.524257 

-0.073113 

0.523527 

-0.072251 

0.522798 

-0.071389 

0.522068 

-0.070527 

0.521339 

-0.069665 

0.520609 

-0.068803 

0.519879 

-0.067941 

0.519150 

-0.069231 

0.517522 

-0.073002 

0.514853 

-0.070245 

0.513497 

-0.067489 

0.512141 

-0.064732 

0.510786 

-0.061976 

0.509430 

-0.059219 

0.508074 

-0.056463 

0.506718 

-0.053706 

0.505363 

-0.050950 

0.504007 

-0.048194 

0.502651 

-0.045437 

0.501295 

-0.042681 

0.499939 

-0.039924 

0.498584 

-0.035411 

0.498243 

-0.030112 

0.498357 

-0.024813 

0.498472 

-0.019514 

0.498586 

-0.017877 

0.499150 

-0.019200 

0.500078 

-0.020523 

0.501006 

-0.021846 

0.501935 

-0.023169 

0.502863 

-0.024492 

0.503791 

-0.025815 

0.504719 

-0.027138 

0.505647 

-0.028461 

0.506575 

-0.026845 

0.507087 

-0.024664 

0.507518 

-0.022483 

0.507950 

-0.020302 

0.508382 

-0.018121 

0.508814 

-0.015940 

0.509246 

-0.013759 

0.509677 

-0.011155 

0.509870 

-0.007495 

0.509467 

-0.003835 

0.509065 

-0.002872 

0.508905 

-0.005687 

0.509086 

-0.008502 

0.509268 

-0.011316 

0.509449 

-0.014131 

0.509630 

-0.016946 

0.509811 

-0.019761 

0.509993 

-0.022575 

0.510174 

-0.025390 

0.510355 

-0.027153 

0.510056 

-0.028465 

0.509551 

-0.029776 

0.509046 

-0.030125 

0.508034 

-0.028444 

0.505953 

-0.026763 

0.503871 

-0.025082 

0.501790 

-0.023401 

0.499709 

-0.021720 

0.497627 

-0.020039 

0.495546 

-0.018366 

0.493471 

-0.018338 

0.492713 

-0.018310 

0.491955 

-0.012521 

0.488779 

-0.007528 

0.486826 

-0.003981 

0.486354 

-0.000434 

0.485883 

0.003112 

0.485411 

0.006659 

0.484939 

0.010205 

0.484467 

0.013752 

0.483996 

0.018188 

0.484059 

0.025236 

0.485693 

0.032283 

0.487327 

0.039330 

0.488961 

0.046378 

0.490594 

0.053425 

0.492228 

0.058037 

0.493241 

0.061599 

0.493987 

0.065160 
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Appendix  E 

Sample  Wave  Parameter  File 


STWAVE  model  result  parameters  Hmo,  TP,  and  mean  wave  direction,  0m, 
are  written  to  file  selhts.out  for  each  selected  grid  cell  identified  in  the  STWAVE 
input  model  parameters  file  (or  SMS  monitoring  station).  A  description  of 
selhts.out  is  given  in  Chapter  4,  “Selected  Wave  Parameters.”  An  example  of  the 
selected  wave  parameter  file  for  Yaquina  Bay  Entrance  at  Newport,  OR 
(Chapter  5,  “Identify  Monitoring  Stations”)  is  given  as  follows: 


Event  ID 

I 

J 

Hmo 

!e 

Ctm 

Breaking 

200700 

110 

157 

8.52 

1.0 

18. 

0 

200700 

131 

117 

7.13 

1.0 

20. 

0 

200700 

146 

86 

4.85 

1.0 

17. 

0 

175500 

110 

157 

4.39 

1.0 

18. 

0 

175500 

131 

117 

4.98 

1.0 

19. 

0 

175500 

146 

86 

4.63 

1.0 

16. 

0 

Appendix  E:  Sample  Wave  Parameter  File 
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Appendix  F 

Sample  Wave  Breaking  File 


Because  the  STWAVE  model  parameter  IBREAK  was  set  to  1,  wave- 
breaking  indices  were  written  to  the  pro  j  ect .  brk  file.  A  description  of  the 
file  is  given  in  Chapter  3,  “Model  Parameters.”  An  example  of  the  first  three 
rows  (J  =  NJ)  of  the  breaking  file  for  Yaquina  Bay  Entrance  at  Newport,  OR 
(Chapter  5,  “Model  Control”)  is  given  as  follows  (0 =non-breaking,  1= breaking ): 
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236 

25.1 

00000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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Appendix  F:  Sample  Wave  Breaking  File 
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Appendix  G 

Sample  Radiation  Stress 
Gradient  File 


Radiation  stress  gradients  were  calculated  and  written  to  the  STWAVE 
output  file,  radstress.  A  description  of  the  file  is  given  in  Chapter  3,  “Model 
Parameters.”  An  example  of  the  radiation  stress  file  for  Yaquina  Bay  Entrance  at 
Newport,  OR  (Chapter  5,  “Model  Control”)  is  given  as  follows: 


236  25.00000 


188 

200700 

6.6950223E-03 
-3 . 3878328E-04 
-7.9356004E-03 
1 . 2184905E-03 
2 . 1935655E-03 
1.2731552E-03 
-9.6608344E-03 
-2 . 0846175E-03 
-5 . 7775499E-03 
-6 . 7550660E-04 
-4.5883558E-03 
-6 . 8725593E-04 
3 . 7447740E-03 
1 . 6133499E-03 
-1 . 1147385E-02 
-2 . 0644760E-03 
-1 . 91 64963E-02 
-3 . 6441803E-04 
-1 . 3087196E-02 
2 . 9396438E-03 
-8 . 0078887E-03 
1 . 0069275E-03 
-6 . 3802721E-03 
8 . 8541035E-04 
-8 . 1194881E-03 
-1 . 3139534E-03 
-1 . 4443645E-02 
-5 . 7125473E-03 
-1 . 7505303E-02 
-8 . 3423238E-03 


-3 . 2863615E-04 
-3 . 9071273E-03 
2 . 5077819E-04 
3 . 54  99572E-03 
2 . 9315948E-03 
-8 . 8462830E-03 
-I . 7865753E-03 
-7 . 9442980E-03 
-1 . 6812515E-03 
-6 . 4390944E-03 
7 . 7571871E-04 
3 . 5375594E-03 
4 . 8  939516E-03 
-1 . 9789848E-02 
-1 . 4142228E-03 
-3 . 7387086E-03 
-2 . 0566939E-03 
1 . 2343216E-03 
3 . 6909487E-03 
7 . 5902939E-03 
3 . 0529022E-03 
-4 . 4814  683E-03 
2 . 3565674E-03 
-1 . 4284382E-02 
-1 . 7316820E-04 
-1.4  678364E-02 
-3 . 3028983E-03 
-1.4  620438E-02 
-6 . 1386488E-03 
-1 . 4727631E-02 


1 . 4794  922E-03 
-1.1065673E-03 
-4 . 3723295E-03 
1.8328476E-03 
-3 . 0966951E-03 
3 . 7113187E-04 
-7 . 0697023E-03 
-1 . 9011307E-03 
-1.0883332E-04 
-2.0248417E-04 
-1 . 1 93801 9E-02 
1.0602188E-03 
3 . 0456923E-03 
-2 . 3662567E-04 
-3 . 6725998E-03 
-3 . 2414244E-03 
1 . 0840988E-03 
2 . 4463271E-03 
-1 . 2021370E-02 
3.2691956E-03 
-8 . 8587953E-03 
1 . 4913178E-03 
-1 . 4439583E-02 
6 . 8014144E-04 
-1 . 0609760E-02 
-1 . 8688391E-03 
-2 . 28  68538E-02 
-5 . 445594  9E-03 
-1 . 8068656E-02 
-8 . 4551238E-03 


1 . 9378655E-05 
-9 . 3047712E-03 
6 . 8103790E-04 
7 . 7973939E-03 
1 . 9795229E-03 
-9 . 5792385E-03 
-2 . 2388841E-03 
-1 . 07  94334E-02 
-8 . 8813779E-04 
-8 . 2439808E-03 
1.9935609E-04 
7 . 97  99647E-04 
3 . 3441163E-03 
-1 . 6355362E-02 
-1 . 7412185E-03 
-2 . 1084480E-02 
-1 . 6369247E-03 
-6 . 6511156E-03 
3 . 0922699E-03 
8 . 6608884E-04 
1 . 584014  9E-03 
2.574  6538E-03 
1 . 6724016E-03 
-9 . 8614125E-03 
-6 . 1447138E-04 
-9 . 32874  62E-03 
-4 . 7854804E-03 
-1 . 5631257E-02 
-7 . 2951126E-03 
-1 . 04  61082E-02 


-2.2567750E-04 
-3 . 7914276E-04 
-2 . 3702241E-03 
2 . 7407454E-03 
-6 . 3614273E-03 
-8 . 2778931E-04 
-5 . 2364729E-03 
-1.9996644E-03 
-7 . 3986058E-04 
5 . 8456417E-04 
-5 . 6306077E-03 
5 . 0053787E-03 
-7 . 4355318E-03 
-1 . 1604691E-03 
6 . 1128233E-03 
-3 . 3995439E-03 
7 . 4631879E-03 
3 . 8751219E-03 
-3 . 4133913E-03 
3 ♦ 5405350E-03 
-1 . 2134209E-02 
2 . 0577621E-03 
-1 . 7014274E-02 
2 . 8734206E-04 
-1 . 6241 912E-02 
-2 . 2675325E-03 
-1 . 9670829E-02 
-5 . 3852084E-03 
-1 . 8562546E-02 
-7 . 3479461E-03 
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